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Project Overview 

The primary goal of the HNC experiment was to obtain high resolution composition 
measurements for cosmic ray nuclei in the platnum-lead and actinide region of the periodic 
table. Secondary objectives include studies of selected groups of elements of lower charge. 
These goals were to be realized by orbiting a large area array of dielectric nuclear track 
detectors be in space for several years. In this time we would collect sufficient actinide 
nuclei to determine the nucleosynthetic age of the cosmic radiation and the relative mix 
of r- and s-process elements in the cosmic ray source. 

The detector was to consist of approximately 50 trays assembled in pressurized can- 
isters. Each tray would contain 8 half-stacks (4 stacks total) and an event thermometer 
which would record the temperature of each event at the time of exposure. Each stack 
would contain seven layers of Rodyne-P, CR-39 and Cronar plastic track detectors inter- 
leaved with Copper stripping foils. 

Upon return to earth, the detectors will be removed for analysis. Ultraheavy nuclei 
would have left “tracks” through the detector sheets that would be made visible after 
etching in a hot sodium hydroxide solution. 

The primary responsibility of the University of Michigan under was to produce and 
assemble the plastic track detectors for the HNC experiment project. The tragic loss of 
the Space Shuttle Challenger first caused an indefinite delay and later a cancellation of the 
HNC project. During the final years, it was decided to use some of the remaining HNC 
funds to upgrade the University of Michigan automatic microscope system that would be 
used to measure large numbers of nuclear that would be returned from such a mission. 

This report contains two parts. The first is a log that describes the activities carried 
out to produce the detector stacks for the HNC experiment. The second part describes 
the automatic microscope system that was developed to scan and measure nuclear tracks. 
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Part I 


Fabrication of the HNC detectors — Project Log 
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Laboratory Set-Up 


The following measures were taken to guarantee security and quality control during 
the detector stack fabrication and assembly: 

B anded stores; Room 148A of the Randall Laboratory' of Physics at the University 
of Michigan was designated for pre-flight storage of all stack materials except CR-39 (see 
bonded cold stores below). Bonded stores is secured by an 8ft. chain-link fence padlock. 

Bonded co l d stores! Room 148 contained an 8 f x 8 ; x 8 1 x 7 " Norlake walk-in freezer 
with holding temperature of — 14°F, to guarantee CR-39 detector sheet preservation over 
pre-flight storage periods. CR-39 sheets were packaged and placed in storage immediately 
after inspection. Bonded cold stores is secured by a lock on the freezer door. 

Light filters: \ ellow light filters were installed over all florescent tubes in the sheet 
handling area (Room 148) to block damaging UV radiation. 

C l eanroom; A filtered clean-room within bonded stores (Room 148A) ensured a dust- 
free environment for the stacking operations (see Appendix E for clean-room layout). 
The cleanroom contains the following equipment for use during the receipt inspection, 
serializing and stacking procedures: 

— IBM PCXT, monitor and printer 

To run coordinating stacking program (see below) 

— Sartorious scale and windshield 

To guarantee detector sheet uniformity by weight (see Appendix C for tolerances). 

Tektronix Digital Plotter and diamond scribe 

For sheet serializing (see Appendix H for scribe design). 

— Heat sealer, aged dry air, Marvelseal packaging material 

For packaging completed stacks and sheet ma terials lft alr/SVat^r/ tl ght tight sea l ' ~ 

— Cleanroom operator garb including nylon lint free hats, smocks, gloves, booties and 
particle masks. 

Stacking program: The stacking procedure was coordinated by the HNC serializing 
and stacking program “BLDSTK.BAS” operable on the IMB PCXT. The program leads 
the operator through the procedure in a step-by-step fashion. Program includes: 

— log m of operator, date and time 

— inventory of available sheet materials to guarantee proper set-up 
weighing procedure to discard sheets outside set tolerances (see Appendix C) 

— scribing procedure to serialize sheets based on serializing code (see Appendix A) 

— stacking procedure to stack sheets in proper order (see Appendix B for stacking 
configuration) 

— packaging and labeling procedure using the plotter to print identifying labels 

— documentation including: 



— continuous printout of stacking procedure 

— floppy diskette data files with sheet serial numbers and weights 

CR-39 Detector Sheet Manufacture 

CR-39 detector sheets were manufactured from monomer with greater than 90% diallyl 
diglycol carbonate (DADC) into 10 mil thick sheets which are laser cut to size following 
NASA blueprint #LC-819728. 

Manufacturing Log: 

— Monomer supplied by SFOS under the trade name “CAD Monomer Allymer” 

5-55 gallon drums received at UM 1/4/85 (UM P.O.#L03450, $1318.75 per drum) 

— CAD monomer tested 

— by UM for viscosity 2/4-2/13/85 (see Appendix K) 

— by Noury Chemicals for chemical and physical properties 4/4/85 (see Appendix J) 
“ two one-gallon samples taken from each monomer drum were stored in Room 148A 

— 1% DOP (dioctyl phthalate) and 0.01% N-445 (Nauguard 445 high performance 
amine antioxidant) by weight was added to the intermixed monomer 

— Monomer and additives were returned to drums and shipped to American Acrylics 
and Plastics (AA&P) to be manufactured into CR-39 sheets 5/6/85 

— CR-39 received from AA&P: 

— 320 sheets (UM P.O.$L00921, $60 per sheet) 

manufactured from standard Barberton monomer with 1%D0P, not CAD-M 

— 44 crates received 6/24-4/14/86 (UM P.O.#L07673) 

100 sheets per crate 

— CR-39 laser cut by Laser Machining, Inc.: 

— 120 sheets, 7/3/85, UM P.O.#L19876 

— 220 sheets, 10/24/85, amendment to above P.O.# 

— 14S sheets, 10/28/85, sent to LaRC for Q/T sheets 
~ 220 sheets, 11/8/85, for Calibration sheets 

— 35 sheets, 11/11/85, for Calibration tabs 


ftodyne-P Detector Sheet Manufacture 


Rodyne-P detector sheets were manufactured from Mobay bisphenol-Apolycarbonate 
resin into 5 and 10 mil thicknesses which are die-cut to NASA blueprint #LC-819727 

specifications. Some sheets contain 0.01% by weight N-512 antioxidant and most contain 
OOP. 

Manufacturing log: 

— Rodyne-P supplied by Coburn Plastics Corp: 

- 500 yards on 20 rolls received at UM 1/25/S5 (UM P.0 #L00313) 

— 7 rolls 20 mil, with N-512 

- 5 rolls 20 mil, without N-512 

— 4 rolls 10 mil, with N-512 

— 4 rolls 10 mil, without N-512 

- 13,175 sheets received at UM 12/4/85 (UM P.O.#L31158) 

— 12,000 sheets 10 mil 

- 1 lot (1300-1500) sheets 5 mil 

— 1100 sheets of 10 mil Rodyne-P die cut by Coburn Plastics Corp. for test stacks 
received at UM 2/6/85 

— GOO sheets 10 mil Rodyne-P die cut by UM Physics instrument shop for Q/T stacks, 
received 12/5/85 



Cronar Detector Sheet Manufacture 


Cronar detector sheets were manufactured by DuPont in 4 mil thicknesses as CR42 
polyester graphics art film which is die cut to NASA blueprint #LC-S19730 specifications. 
Manufacturing log: 

— Cronar supplied by Northern Graphics supply 
13,500 sheets (133 boxes) received at UM 11/1/85 

— Cronar sheets die cut for Q/T stacks by UM physics instrument shop 

— 400 cut sheets received 12/5/85 

— additional cut sheets received 4/22/86 

Copper Sheet Manufacture 

Copper sheets were manufactured from type-101 OFHC Certified copper ASTM F6877 
into 10 mil thick sheets which were die-cut to NASA blueprint #LC-819729 specifications. 
Manufacturing log: 

— Copper supplied by Olin Brass 

5900 lbs. (approx. 6775 sheets) received at UM 4/14/86 (UM P.O.#C97746, approx. 
$15,197.75) 

— 1 case copper sheets die cut for Q/T stacks by UM physics instrument shop, received 
4/22/86 

Calibration Sheets and Calibration Tabs 

Calibration sheets are plastic track detectors exposed to relativistic 56 Fe nuclei (for 
CR-39 sheets) or “ 38 U nuclei (for Rodyne-P and Cronar sheets) at the Lawrence Berkeley 
Bevelac. They were to be used to monitor post flight thermal annealing of tracks. 

Calibration tabs were detectors irradiated with the calibration sheets. They do not fly 
with the experiment but remain in storage until experiment completion for comparative 
analysis and identification. 

Calibration packets each containing enough calibration sheets and calibration tabs for 
one canister were packaged in Marvelseal and sent to Bevelac at the University of California 
in Berkeley for irradiation. 

See Appendix B for exact placement of calibration sheets in the detector stacks. 

Qualification and Test Stack Assembly 

Qualification and test (Q/T) stacks were built to test the overall detector and canister 
design for the HNC project. They are assembled identically to proposed detector stacks 
except that unirradiated detector sheet is substituted for each calibration sheet in the stack 
configuration. The completed stacks were sent to the NASA Langley Research Center for 
testing. 



HNC Project Mothballing 


Due to the space shuttle Challenger accident, further work on the HNC project was 
originally postponed until a new flight date could be scheduled. To preserve the detectors 
until such time, all plastic sheet materials were double-bagged in Marvelseal packaging 
material and purged with dry air before being sealed and placed in bonded cold stores. All 
copper sheets remained crated in bonded stores. 
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33 

rejects 

100 

2-20-86 

34 

3-3-86 

100 

3-4-86 

35 

3-31-86 

100 

4-8-86 

36 

4-14-86 

100 

5-1-86 

37 

5-1-86 

100 

5-1-86 

38 

8-8-86 

100 

8-12-86 

39 
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8-12-86 

40 

rejects 

73 

8-12-86 
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cal shts./tabs 

125 


crate 2 

cal shts./tabs 

125 
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Rodyne-P 


Packet No. 

Date Pkgd. 

No, of Sheets 

date in cold stores 

thickness (mil) 

0001 

3/5/87 

125 

3/13/87 

10 

0002 

3/5/87 

150 

3/13/87 

10 

0003 

3/5/87 

125 

3/13/87 

10 

0004 

3/6/87 

125 

3/13/87 

10 

0005 

3/9/87 

125 

3/13/87 

10 

0006 

3/9/87 

125 

3/13/87 

10 

0007 

3/9/87 

125 

3/13/87 

10 

0008 

3/10/87 

125 

3/13/87 

10 

0009 

3/10/87 

125 

3/13/87 

10 

0010 

3/10/87 

125 

3/13/87 

10 

0011 

3/10/87 

125 

3/13/87 

10 

0012 

3/10/87 

125 

3/13/87 

10 

0013 

3/10/87 

125 

3/13/87 

10 

0014 

3/10/87 

125 

3/13/87 

10 

0015 

3/10/87 

125 

3/13/87 

10 

0016 

3/10/87 

125 

3/13/87 

10 

0017 

3/10/87 

200 

3/13/87 

5 

0018 

3/10/87 

200 

3/13/87 

5 

0019 

3/10/87 

200 

3/13/87 

5 

0020 

3/10/87 

200 

3/13/87 

5 

0021 

3/10/87 

170 

3/13/87 

5 

0022 

3/11/87 

125 

3/13/87 

10 

0023 

3/11/87 

125 

3/13/87 

10 

0024 

3/11/87 

125 

3/13/87 

10 

0025 

3/11/87 

125 

3/13/87 

10 

0026 

3/11/87 

125 

3/13/87 

10 

0027 

3/11/87 

125 

3/13/87 

10 

0028 

3/11/87 

125 

3/13/87 

10 

0029 

3/11/87 

125 

3/13/87 

10 

0030 

3/12/87 

125 

3/13/87 

10 

0031 

3/12/87 

125 

3/13/S7 

10 

0032 

3/12/87 

125 

3/13/87 

10 

0033 

3/12/87 

125 

3/13/87 

10 

0034 

3/12/87 

125 

3/13/87 

10 

0035 

3/12/87 

125 

3/13/87 

10 

0036 

3/12/87 

125 

3/13/87 

10 

0037 

3/12/87 

125 

3/13/87 

10 

0038 

3/12/87 

125 

3/13/87 

10 

0039 

3/12/87 

125 

3/13/S7 

10 

0040 

3/12/87 | 

150 

3/13/87 

10 
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3 / 12/8 

3 / 12/8 

3 / 13/8 

3 / 13/8 

3 / 13/8 

3 / 13/8 

3 / 13/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 19/8 

3 / 20/8 

3 / 20/8 

3 / 20/8 

3 / 20/8 

3 / 20/8 

3 / 20/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 23/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 

3 / 24/8 


125 

3 / 13/87 

10 

150 

3 / 13/87 

10 

150 

3 / 13/87 

10 

125 

3 / 13/87 

10 

150 

3 / 13/87 

10 

150 

3 / 13/87 

10 

150 

3 / 13/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 

125 

3 / 25/87 

10 
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0084 


3/24/87 

0085 


3/24/87 

0086 


3/24/87 

0087 


3/24/87 

0088 


3/24/87 

Misc. A 


3/25/87 

Misc. B 


3/25/87 

Misc. C 


3/25/87 

Misc. D 


3/25/87 

Misc. E 


3/25/87 

Misc. F 


3/25/87 

Misc. G 


3/25/87 

Misc. H 


3/25/87 

Misc. O 

L 

3/25/87 


M 

3/25/87 


N 

3/25/87 


125 

125 

125 

125 

125 

150 

200 (QT shts) 
approx. 150 ( thick. vax.) 
150 

approx. 100 (thick.var) 
150 
150 
150 

100 (some cut,?) 
approx. 100 (some ?) 
approx. 100 (some ?) 


3/25/87 

10 

3/25/87 

10 

3/25/87 

10 

3/25/87 

10 

3/25/87 

10 

3/26/87 

10 

3/26/87 

5 

3/26/87 

10 

3/26/87 

10 

3/26/87 

10 

3/26/87 

10 

3/26/87 

10 

3/26/87 

10 

3/26/87 

5 

3/26/87 

10 

3/26/87 

10 
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Cronar 


Packet No. 

Date Pkgd. 

No. of Sheets 

Date Into Cold Stores 

0001 

2/9/87 

300 

3/2/87 

0002 

2/9/87 

300 

3/2/87 

0003 

2/9/87 

300 

3/2/87 

0004 

2/9/87 

300 

3/2/87 

0005 

2/9/87 

300 

3/2/87 

0006 

2/9/87 

300 

3/2/87 

0007 

2/9/87 

300 

3/2/87 

0008 

2/11/87 

300 

3/2/87 

0009 

2/11/87 

300 

3/2/87 

0010 

2/11/87 

300 

3/2/87 

0011 

2/11/87 

300 

3/2/87 

0012 

2/11/87 

300 

3/2/87 

0013 

2/12/87 

300 

3/2/87 

0014 

2/12/87 

300 

3/2/87 

0015 

2/12/87 

300 

3/2/87 

0016 

2/27/87 

300 

3/2/87 

0017 

2/27/87 

300 

3/2/87 

0018 

2/27/87 

300 

3/2/87 

0019 

2/27/87 

300 

3/2/87 

0020 

2/27/87 

300 

3/2/87 

0021 

2/27/87 

300 

3/2/87 

0022 

2/27/87 

300 

3/2/87 

0023 

2/27/87 

300 

3/2/87 

0024 

3/2/87 

300 

3/9/87 

0025 

3/2/87 

300 

3/9/87 

0026 

3/2/87 

300 

3/9/87 

0027 

3/2/87 

300 

3/9/87 

0028 

3/3/87 

300 

3/9/87 

0029 

3/3/87 

300 

3/9/87 

0030 

3/3/87 

300 

3/9/87 

0031 

3/3/87 

300 

3/9/87 

0032 

3/3/87 

300 

3/9/87 

0033 

3/3/87 

300 

3/9/87 

0034 

3/3/87 

300 

3/9/87 

0035 

3/3/87 

300 

3/9/87 

0036 

3/3/87 

300 

3/9/87 

0037 

3/3/87 

300 

3/9/87 

Misc. I 

3/25/87 

300 

3/26/87 

Misc. J 

3/25/87 

300 

3/26/87 

Misc. I< 

3/25/87 

300 

3/26/87 

Cal. Sheets 

3/3/87 

389 (1200-?bad) 

3/9/87 
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Appendix A 

Sheet Serializing Code 


Calibration sheets are scribed with an 11-digit serial number. All other sheets are 
scribed with a 9-digit number. The codes are as follows: 


Cal-sheet identification 

CS 

(cal-sheets only) 


Canister number designation 

01- 

Stack-quad location identification 

A-D 

Canister side location 

1-2 

Layer numbers 

1-7 

Layer sheet numbers 

01-12 

Stack sheet numbers 

01-72 


Example: 05C 160933 

The above example identifies a hypothetical sheet as located in: 

1) the fifth (5th) CANISTER, in 

2) Quad “C”, on 

3) CANISTER SIDE “One” (1), in 

4) LAYER NUMBER “Six” (6), being 

5) LAYER SHEET NUMBER, “Nine” (09), and is 

6) STACK SHEET NUMBER, “Thirty- Three” (33) 


H 



Appendix B 

HNC Stack Configuration 


No. 

Item 

No. 

Item 

No. 

Item 

No. 

01 

R 

12 

CR 

23 

CU 

34 

02 

R+ 

13 

CU 

24 

C 

35 

03 

CU 

14 

c 

25 

C 

36 

04 

C+ 

15 

c 

26 

C 

37 

05 

c 

16 

c 

27 

C 

38 

06 

c 

17 

c 

28 

CU 

39 

07 

c 

18 

CU 

29 

R 

40 

08 

CU 

19 

R 

30 

R 

41 

09 

R 

20 

R+ 

31 

R 

42 

10 

R+ 

21 

R 

32 

CR 

43 

11 

R 

22 

CR 

33 

CU 

44 


Item 

No. 

Item 

No. 

Item 

No. 

Item 

C+ 

45 

C 

56 

c 

67 

C 

c 

46 

C 

57 

c 

68 

CU 

c 

47 

C 

58 

CU 

69 

R 

c 

48 

CU 

59 

R+ 

70 

R 

CU 

49 

R 

60 

R 

71 

R 

R+ 

50 

R 

61 

R 

72 

CR 

R 

51 

R 

62 

CR 



R 

52 

CR 

63 

CU 



CR 

53 

CU 

64 

c+ 



CU 

54 

C 

65 

c 



c 

55 

C 

66 

c 




CU = copper 
R = Rodyn-P 
C = Cronar 
CR = CR-39 
-I- = skeleton etch 
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Appendix C 

Sheet Weight Tolerances 


program ^BLDSTiTbAS” , ° ler ‘“ CeS " ere ^“rporated into the senalizing and stacking 

ficition hv th K '? AS i ^ detector sheet and proper material identi- 

fication by the operator dunng the stack assembly procedure. 


Material 

10 mil Rodyne-P 
4 mil Cronar 
10 mil CR-39 
10 mil copper 


Iblerance max (in grams) min 


59 

32 

70 

376 


45 

21 

50 

368 
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Appendix D 

HNC Correspondence References 

AKZO-CHEMIE America 
Noury Chemicals 
8401 W. 47th Street 

Attn: Bob, Maxgevich, Technical Development Manager 

American Acrylics & Plastics 
25 Charles Street 
Stratford, CN 06497 
Attn: George Barry, President 

Cobum Plastics Corp. 

1650-T Corporate Rd. W 
Lakewood, NJ 08701 

Attn: Paul Riske, Manager Industrial Films 

Laser Machining, Inc. 

500 Laser Drive 

P.O. Box 219 

Somerset, WI 54025 

Attn: Ted Bauer, Robert Bosworth 

Northern Graphics Supply 
2216 Midland Road 
Saginaw, MI 48603 

Olin Brass 

19 Atterbury Blvd., Suite 11-B 
Hudson, OH 44236 

O’Neal, Bob 
Project Manager 
Mail Stop 258 

NASA Langley Research Center 
Hampton, VA 23665 

Price, Professor Buford 
Department of Physics 
University of California, Berkeley 
Berkeley, CA 94720 

Sola Optical 
3600 Lakeville Hwv. 

Petaluma, CA 94952 
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Appendix H 

Diamond Scribe 


Appendix I 

Temporary Stacking Hardware 


The following hardware was manufactured at the University of Michigan to be used 
during the stacking operation for stabilizing purposes. All temporary stacking hardware 
was removed from the stacking fixtures upon completion of the stacking operation and 
before installation of the mounting hardware. 


A) Aluminum Bushing 



0. D. 0.435 in. 

1. D. 0.325 in. 
quantity: 15 


B) Stabilizing Bolt 



i\n. 
♦ . 

1 

1 

t 




O.D. 0.30 in. 
steel 

5/16 - N.C. 
quantity: 20 


Installation of the Temporary Stacking Hardware (cross-section) 



ORlG-?r-’.U PAGE IS 
OF POOR QUALITY 


Appendix J 


Analysis of CAD-N Monomer 
Akzo Chemie America. Noury Chemicals 
Ref: ARCS-39 

Project Number 06S4075039 

Eight samples of CAD-N Monomer were received from the University of Michigan for 
analysis. These samples were analyzed by wet methods, physical analysis. HPLC, and gas 
chromatography. The results of the wet and physical analyses are shown in Table 1. The 
viscosity was done by the falling ball method. Karl Fischer titration was used to measure 
water. The pH was measured on a 10% solution in 50/50 water-ethanol. Inorganic and 
organic hydrolyzable chlorine was measured using the Noury procedure. 


Table 1 

Wet and Physical Analysis of CAD-N Monomer Samples 



i 

2 

3 

SAMPLE NO. 
4 5 

6 

7 

8 

Viscosity (mpa • s )* 

IS. 5 

18.7 

18.5 

IS. 4 

18.5 

17.9 

18.5 

19.2 

Water (ppm) 

797 

791 

733 

677 

749 

1186 

541 

477 

Color (APHA) 

5 

5 

5 

5 

5 

5 

5 

5 

pH 

5.2 

5.7 

5.5 

5.5 

5.6 

4.8 

4.6 

4.7 

Refractive Index 

1.4530 

1.4522 

1.4522 

1.4521 

1.4520 

1.4525 

1.4517 

1.4527 

Inorganic (ppm) 
Organic 
Hydrolyzable 
Chlorine 

2.5 

2.3 

2.4 

4.2 

4.2 

IS 

ll.S 

9 


+ Dynamic viscosity millipascals x seconds 
Sample Code: 

1-5 = Allymer drums no. 15 

0 = misc. 

7 = #3178-5-16 

S = #5239 #5 
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High performance liquid chromatography (HP LC) was used to determine the monomer 
• hmer and trimer content. The allyl diglycol carbonate was also measured. It was also 
onserved that all samples contained some polymer. The HPLC results are shown in Table 


Table 2 


HPLC Analysis of CAD-N Monomer Samples 



1 

2 

3 

SAMPLE NO. 

4 5 6 

7 

8 

n = 1 

92.4 

92.1 

86.6 

93.6 

8S.0 

93.0 

86.8 

94.1 

n = 2 

4.5 

4.3 

3.8 

4.5 

4.1 

3.9 

5.3 

4.9 

n = 3 

0.6 

0.5 

0.4 

0.3 

0.3 

0.8 

0.8 

0.4 

Allyl Diglvcol 

1.6 

1.6 

1.4 

1.6 

1.3 

1.5 

0.9 

0.5 

Carbonate 










GaS chromat °g ra phy was used to measure the acrolein (ppm), allyl alcohol (ppm) and 
diallyl carbonate (%) content. The results are shown in Table 3. 

Table 3 


Gas Chromatographic Analysis of CAD-N- Monomer 



1 

2 

SAMPLE NO. 

3 4 5 6 

7 

8 

Acrolein (ppm) 

68 

40 

1 7S 73 

49 256 

127 

147 

Allyl Alcohol 
(ppm) 

58 

29 

32 43 

36 88 

73 

94 

Diallyl Carbonate 
(%) 

0.25 

0.21 

0.25 0.21 

0.22 0.31 

0.31 

0.33 





Appendix K 

Viscosity Measurements of CAD - Monomer 


The viscosity of the monomer in each of the five drums received 1/4/85 (UM P.0. 
#L03450) was measured at the University of Michigan prior to the monomer use in the 
manufacture of CR-39 sheets. 

The following data was obtained using the '‘Falling Ball Viscometer Method” (vis- 
cometer accuracy: 0,2 - 1.0%): 


Sample No. 

Centipoise 

Avg. Centipoise 

Denistv 

1 

15.0059 



1 

14.9531 

14.9669 

1.1406 

1 

14.9418 



2 

14.8605 



2 

15.3168 

15.0151 

1.1407 

2 

14.8651 



3 

14.92S8 



3 

15.0082 

15.0325 

1.1408 

3 

15.1107 



4 

14.8005 



4 

14.9124 

14.9014 

1.1411 

4 

14.9914 



5 

14.9731 



5 

15.0597 

15.0532 

1.1412 

5 

15.1268 





The following data was obtained from density measurements used to check the 
t the viscosity measurements: 


Sample No. 

Mass (g) 

Volume (cc) 

1 

114.06 

100 

2 

114.07 

100 

3 

114.08 

100 

4 

114.11 

100 

5 

114.12 

100 

avg. density = 1.14088 


accuracy 




Part II 

Description of the University of Michigan 
Automatic Nuclear Track Measuring Microscope 


1. Introduction 

[“‘fl 18 8eCt ‘°" ' he i ardware s °f‘ware f°r the track measuring microscope svstem 
| 16 5 . eSCn e : Tbe on S mal S°al of the group was to develop a svstem capable of 

ZlZf meaSUrmg 1 milli ° n COSmiC ' ray CVentS automatically, for u sTon a 

. .ety of expenments requiring large collection areas and precise charge resolution While 

erisZr r ^ "“'“‘'“a ° f CR ’ 39 Pl3StlC " Wrfl >“»"• th ' 'Nations imposed bv 
clear tracksZ pr ° Cess,n S ^jechnoiogy make accurate identification and measurement of nu- 

Ub over th "7“ t,me a f °™ da We task. The system that has evolved in our 

b over the past four years represents a compromise between speed and accurate identi 
hcation and measurement suitable for our purposes. The time required for meZurement 

our Zd M t at ‘ nd “f ted above - Whlle farther improvements in hardware could improve 
the “cent L SVh? y ’ a ? Va ” CeS in t‘ he devel ° pment ot Processing systems since 

OtovZ the <■ ^ PrOJeCt (c,rca 1985) hold the P^P" 1 of dramatically im- 

vpl ltm i Y Cq T?a for , analysis “ d software development for experiments of this 

find application 5 in < hUure^rstmns^. eC ^ n ' qUeS and “ d ” " 

deta^n^H f °! CR - 39 d ^ be described briefly here and in more 
n section 4^ This procedure begins by determining a coordinate svstem on the 
^ eet of interest, which will serve as a reference coordinate system for all subsequent track 

ooZ P n ™ c rrd Thi t s c , oordina : e r em is “ d ty tw ° p~-. *. zsz 

oT? at th : poin ; ° f interest - a « *■ zz Zch 

defines the ong n (0.0). In practice, the origin of the coordinate svstem is taken to be the 

“ ° a " *“* * -rmally incident to the surface of the detector The 

potion of a second track is then used to define a reference coordinate svstem /= 0) 

no t on f ,b T f T ‘ he Sta * e ' and ‘ he " repla “ d ' “f station introduced in the 
posit, on of the sheet may be removed by locating the reference tracks and comparing^ 

’T COOrdm *‘« tQ those of ‘he reference locations of these same tracks. In this wa a 
subsequent scans of the sheet may be referenced to identical coordinate systems 

Next each side of each sheet is scanned at low power (typically 3 or 10X) «n<l 
. he localioiis of all candidate tracks are idenlified and their c^Ldniatesthnm ,o disk 

the baT Van TT arame,CrS Whkh have been ralculated for these tracks. In to eliminate 

from T n t SU “ e eCtS “ d non 'Penetrati n g tracks from low energy particles 
from consideration as candidate events, the opposite side of this surface is afto scTimed 

Zh P T,b * e '' entS aK idenlified by roinc, deuce between these two surfaces after tte 
paths of the events ,n the plastic are reconstructed. 

The results of this scan are then used to make final measurements of track parameters 

‘ th P< 7”. " 7° tllCSe mcasurements <t«n ''d made more accurately improve the 
charge resolution that can be obtained from the plastic. ' 

These procedures have been used on data obtained from accelerator experiments and 
also „„ balloon flight data. The performance of the system is discussed in section o 
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2. Microscope Hardware 


A schematic diagram of the automated measuring system is shown in figure 1. The 
microscope consists of a Technical Instrument Measuremaster II stage with 5 Leitz objec- 
tives of 3.2X, 10X, 20X, 32X and SOX magnification. The microscope is equipped with a 
Hamamatsu C1000 Vidicon camera. The output of the camera is directed to a VICOM 
model 1850 image processing system, where the image is digitized and subsequent analysis 
is performed. Mass storage is provided by a CDC 450 Megabyte hard disk, where the 
iesults of the image analysis are collected. Data files can be transferred to and from the 
\ ICOM by means of a DMA link with a DEC micro VAX. where the results can be further 
analyzed and plotted without using CPU time dedicated to image analysis. 

The stage of the microscope system has a 6X6 "travel on 1mm pitch lead screws and 
is controlled by two 25,000 step/revolution Compumotor stepping motors and Compumo- 
tor 2100 series indexers. The position of the stage is monitored by a Heidenhain length 
measuring system, giving absolute position information accurate to ±l//m over its entire 
length of travel. The focus of the microscope is monitored by a Leitz Laser Autofocus 
system. 

The VICOM image processor is built around a Motorola VM03-4 microcomputer, and 
includes a fast point processor and array processor in a pipeline architecture. This allows 
\arious image operations, such as region of interest processing, histogramming and 3X3 
pixel convolution, to be accomplished at the frame rate (approximately 1/30 second). 
These capabilities have proved to be extremely useful in the design of the image processing 
algorithms required to smooth image data and extract edges from potential tracks, as 
described below. The images are digitized by a 8 bit /pixel analog to digital converter. 
The digitized images are written to a 512X512X2 byte segment of memory dedicated to 
image storage. Three additional segments of memory are available to store the result of 
image processing operations. The VICOM also has some graphics capabilities, with 4 bits 
per word in image memory reserved for graphics display. 

The microcomputer also controls the stage and autofocus via an RS-232 link with the 
Compumotor indexers. The motion of the stage is initiated by serial command strings 
sent by the microcomputer. The indexers contain their own GS02 microprocessor, so the 
operation of controlling and monitoring the position of the stage is divorced from the 
\ ICOM microprocessor itself, and real-time analysis on digitized images can proceed while 
stage movement is in progress. 

The Leitz Laser Autofocus system can be controlled via an operator keypad or by 
means of programmable bits on the Compumotor indexers which can be accessed from the 
RS-232 link. In practice, both methods of control are used during the plastic analysis, 
as discussed below. The focus axis is also equipped with a Heidenhain length measuring 
system and special electronics to give .1 /ira position resolution over the entire length of 
tia\el. Optionally, the position of the stage can be controlled by a Compumotor stepping 
motor to this same accuracy. 
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During analysis, the plastic is mounted on a specially constructed aluminum stage 
■ quipped with a vacuum line that holds the sheet in place during analysis and keeps the 
sheet flat to within .001”. This design restricts us to the use of reflected light for viewing 
samples. When viewing samples using a low power, long focal length objective, a reflective 
paper backing is used to increase image contrast. This also improves the performance of 
the laser autofocus by reducing reflection from the laser off of the aluminum stage. 

The entire microscope apparatus is mounted on a large marble block and isolated from 
vibrations by a specially designed microscope stand with sand-filled legs. 
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3. Analysis Software 


The Motorola VM03-4 microcomputer in the model 1850 VICOM image processing 
system runs under the VERSAdos operating system and includes a resident PASCAL 
compiler and Motorola assembler for 68000 based systems. VICOM provides a library of 
general-purpose image processing routines and software necessary to program the special 
purpose hardware. 

One of the original goals in developing the microscope software was to provide an 
environment for quickly setting up and measuring sheets by an essentially unskilled worker 
with a minimal amount of training. Another goal was to make the program entirely 
modular to save time in development of the routines required for the plastic analysis. To 
these ends we have designed the system around a central menu program, which provides 
an easily learned user interface as well as the capacity to load and run the various analysis 
programs and monitor their status. This structure allows the development of the different 
analysis programs to occur simultaneously, and avoids the problems inherent in developing 
all of the analysis routines in a single monolithic structure. While menu system remains 
resident in memory, the data analysis programs can be swapped in and out of memory 
under control of the menu and the limited memory space on the VICOM is very efficiently 
used. 

The menu user interface is mouse driven and can be operated on a Tektronix 4014 
compatible terminal. The VERSAdos operating system is built around a real-time multi- 
tasking kernel, and offers wide variety of sophisticated intertask communication facilities 
suitable for our purposes. The menu also contains a shared segment of memory containing 
the data structures necessary for passing information between analysis programs. Figure 2 
shows a schematic diagram of the menu system and associated analysis tasks. 

The heart of our analysis system are the image processing routines. The accuracy 
with which these identify and measure tracks largely determines the resolution which 
the experiment will have, and the speed at which these routines execute is the primary 
constraint on total analysis time. Our philosophy throughout has been to dedicate the 
VICOM entirely to low-level image processing tasks, while data reduction and further 
analysis done on the microVAX. The problems that this has introduced into our analysis 
procedure are discussed below. The VICOM image processor contains a large library of 
PASCAL-callable subroutines and a resident PASCAL compiler, and thus we have chosen 
to program the analysis programs in PASCAL, and to use assembly language for portions 
of the programs where speed is essential. 

As discussed previously, the microscope image is produced in reflected light. Tracks 
produced with various angles of incidence appear as ellipses with various internal features 
which generally reflect considerably less light than the horizontal plastic surface. The goal 
of our routines is to convert these images to a list of digital information which can be 
processed to give us the dimensions and orientation of the ellipse which best approximates 
the perimeter of each track. Since all of the necessary information about each track (major 
and minor ellipse axes, and ellipse orientation) can be derived from its perimeter points 



P ° in !: T titUte a Smaller and m ° re mana S eable ■* of data than the 
1 ixe s that make up the ellipses area, we have concentrated on developing edge enhance 

aU 7 th f P~r to be identified and scented from the 
gr und of the image. An alternative scheme, which involves the segmentation of the 
Pses by using an optimal threshold selection (on the basis of a histogram of the entire 
• mage) was a empted m the early phases of the development of our system. This method 
'A as not considered practical because a) the tracks are not uniformly dark, and therefore a 

r USCd t0 SCleCt the track fr ° m ba <*ground will not faithfully 

wh b V he u the track ’ and b) the reSUlt ° f the thresh °Wing is a list of pixels 
ich make up the interior points of the track ellipse, and, as mentioned above this is a 

arger set of pixe s than are needed to evaluate the ellipse parameters. Furthermore algo- 

n W used to stdect the edge points front this se, of pixels are computationally expend 

1 6 ’ ‘ heSe alg0nthms ca " ba implemented as spatial operations on a single image. 

There is a considerable literature on the subject of detecting edges in machine vision 

d™ 7auT„ ^h f 5i Pdi “1 MiJah 1982; Leraldi 1983 > “ d • i-ge class of these 

nu hes In ™ enhancement/threshold category. In these edge detectors, disconti- 

the rnh d T**' nbU ‘ e are enhanced or accentuated by some spatial operator. If 
In cdce d T‘ m “ ty sufficiently ‘"Se, that is, greater than some threshold level 
Lrah e f ^ p ^"‘( figure 3). This type of edge detection scheme particularly is 

f ew of th” 1 POm ‘ ° f I f t haidW “ e -d from the pit 

of v e„ of the software required for feature-list extraction and subsequent analysis As 

conv ? n . ed PreVKWsly ’ our ,ma S e Processing facilities include the capabilities for point and 

“"a tsuuTr T tQ e Pe Tr d qUkkly - “ d lh ' . t hat must belalld 

An obZ V t, enhancement/threshold operation ellipse parameters is relatively small. 

Ksclr , “ he SPa,ial ° peral ° r iS a « ««<* derivative operation 

fhe edge ” ° ,mage ,ntenSHy aPPear “ " Peak ” grey ICVd a ‘ ‘ he l0Cati " n of 

HiIdretlW198nl' Ve ed * es in ima 8“ has been proposed by Marr and 

eldlion of he T' “ d POggi ° (1979) “ * p “ sibk mechanism for the 

raw primal sketch in early visual systems. This scheme which is 

presumably implemented neurally in the visual system, uses a two-dimensional extension 

of t b^ D l" 0 «n aS ri L ? g “ . S the ° re " 1 (L °«a" 1!)77 )- I‘ states that the se«> crossings 
a band-pass filtered function can be used to determine the original function uo to a 

mu tiphcat^e constant, provided the function satisfies some conditions for being '‘‘well- 
behaved . In this method, which is discussed more fully in appendix A, theLI is 

Z° The result fflh' Wl " Ch 8 Gaussian smoothing operation with a Laplacian 

nass filter Th 1 f hlS sequence of opcmtions is an approximation of a one octave band- 
pass filter The zero crossings are then detected and identified with edges. While the 

heortd, cal motivations for this method of edge detection are not entirely sound it I 

which t '’T Se " " 0t rely ° n a threshold is therefore applicable to images 

r„dth d CO " lraSt '™ ali0n ' A|S °' the La P ,acia " « 1( ct is spherically symmetric 

and the edges which are extracted do not show the directional dependence apparent in some 
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edge detection schemes that involve directionaly dependent spatial derivative operators. 

The edge detection algorithm that we have implemented is an enhancement/threshold 
operation which incorporates the ideas presented by Marr et. al. Although we have not 
implemented the zero crossing detection that this procedure suggests, the thresholding 
scheme that we use for determining the gray level peaks in the second order derivative 
of the smoothed image is adequate for the purposes of interest, while remaining easy to 
implement in software. We have not used convolutions which are identical to those of 
Marr et. al., but instead have chosen a set of 3X3 convolutions which correspond to a 
spatial smoothing filter followed by a Laplacian operation. The necessity of the smoothing 
becomes apparent by studying the results of a Laplacian operation on an image with and 
without smoothing. This smoothing tends to eliminate noise in the image which has a scale 
smaller than the size of the 3X3 pixel mask, such as digitization errors, or very small scale 
imperfections in the plastic surface, while preserving large scale features, such as track 
edges (see appendix A). This reduces the background noise in the thresholded image, and 
simplifies the task of segmenting tracks from this background. 

Our first attempt at using the ideas presented above used the convolution masks shown 
in figure 4, followed by a gray level threshold to produce a binary edge map in image 
memory. The image containing this map was then searched for the occurrence of an “on” 
pixel, and a list of all 8-neighbor continguous pixels was found (using a depth-first search 
algorithm) and stored in a linked list structure, which was then used to compute the ellipse 
parameters. This algorithm worked well if the track edge was well defined, but failed for 
certain commonly occurring track imperfections. This tended to break the continuity of 
the track, making it difficult to identify contiguous segments of pixels with single tracks. 

The next step in the evolution of the edge detection algorithms was the implementation 
of a threshold/tracking detector (appendix B). The program is conceptually very simple. 
After the smoothing/Laplacian operations have been performed on an image, the pixel 
rows are scanned for the occurrence of a single pixel above a certain gray level. After a 
pixel has been found, its position is recorded on a stack, and a list of adjacent pixels are 
scanned for a value above some second threshold. If a pixel is found that satisfies these 
conditions this pixel is stored a linked list of track candidate pixels. If more than one 
pixel is found above threshold, only the one with the highest gray level is accepted. This 
procedure then continues at the location of the pixel most recently found. The search ends 
when a) the location of the next pixel is equal to the start position of the search, in which 
case the linked list is terminated, and a track is assumed to be found, or b) the search 
reaches a “dead end”, that is, the start position has not been found, and no adjacent pixels 
satisfy the threshold condition for the search to continue. In this latter case, the linked 
list of pixels for the candidate track is deleted from memory. Thus each track has a list 
of pixels which forms a continuous ring, which helps to remove from consideration edges 
formed by objects which are not tracks (an example is a surface scratch in the plastic 
which appears as a long straight edge). A track candidate is then tested to make sure that 
the number of pixels found lies above some certain value; if the track does not pass this 
cut the linked list is also deleted from memory. Control is then returned to the section of 
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the program that scans lm age memory for pixels above threshold. In practice, to speed up 
ns operation, only select rows of the image are scanned; the spacing between these rows 
lire chosen such that they are smaller than the size of the objects of interest. 

In the process of performing this search the x and y coordinates of the pixels are 
summed to provide first-order moments for calculating the centroid of the track. All of 

speed. 1006 11X6 ^ WrittCn m °P timize d assembly language to provide the fastest possible 

After the entire image has been scanned in the manner described above the linked lists 

of candidate tracks is passed to a PASCAL routine which does the subsequent processing 

The primary reason for using the linked list as opposed to an array structure is that the 

i ed list is more flexible when dealing with pixel lists of varying size. For example if an 

array structure were used, the arrays would have to be made large enough to accommodate 

tracks viewed at high power, while at low power, most of this memory would be wasted 

Also, many of these arrays would have to be used to store the results at low power, when 

a large number of tracks may appear in an image. When measuring at high power, on 

e o er and, in general only one track is of interest to the analysis. The linked list 

PA^A T “*! lmplemen ‘ in ^en>bly language; this type of data structure appears 

fron, PA en “ re Stn,CtUre can be instructed in assembly language and accessed 

irom i AoLAL, as necessary. 

The parameters of interest in the surface analysis of the tracks are the major and 
minor axes of the ellipse forming the track mouth, and the orientation of this ellipse 
These quantities, along with the etch rate of the bulk surface of the plastic iv suffice to 
determine the charge of the particle, given that the particle is relativLic 

These parameters are calculated using a least-squares fit to the pixel list. The method 
is given in appendix C. The parameters are calculated from the ellipse moments: 


M 


mn 


E 


r"V 


where * and y are the x and y coordinates of a pixel and the sum runs over each pixel in 
a rac he derivation of the ellipse parameters from the moments for the most general 
case is quite complicated and requires 15 moments to be calculated. In practice ive use an 
approximation to this general solution, which requires only 10 moments for all the param- 

r<T J °”" d mm ° r “?’ orienUtion “Sk) “> b “ calculated. The relative complexity 
of the 10 and 15 moment solutions can be inferred from appendix C, where the solutions to 

gra^ REDUCE ST T. '™* ^ “*“* ““ the algebraic manipulation Pro- 

Sa 3Xw ? ,r a S,mphc,ty of the 10 moment s°lP«°n (requiring the solution 

i e yS em , ° f mear lnhom °«cneous equations, compared with a 5 X 5 system for the 

pmdXforX fir.^r' b! ' ‘ he feCt ^ hat the Cen,r ° id ° f ‘ he ellip5e must be d«ennined 
usTnHh A/ 1 ,°, MCUrate ” this nPPtoximation. Since the centroid is determined 
it e S t 1 wf T monients of the elbpse - the determination of the centroid makes 

h , ,1 ,v * f ^mentation algorithm find the entire edge of the ellipse, and 
that this edge have uniform thickness. A portion of the edge of an ellipse, for example 
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would give a bad fit in the 10 moment case, since the centroid of the full ellipse would not 
be accurately determined by the moments derived from the portion of the ellipse found. 
The threshold/ tracking algorithm that we have implemented satisfies both of these con- 
straints: the edge is of single pixel width, and the completeness of the edge is assured by 
the requirement that the last pixel in the search be the pixel from which the search was 
started. 

In order to increase the speed of the analysis, none of the parameter calculations are 
done on the VICOM at low power. Instead, the moments are summed and written directly 
to disk, along with the sheet x and y coordinates of the track. Each image is also assigned 
a number and each track candidate within the image assigned a unique ‘event number’. 
These moment files are then transferred to the micro VAX, where ellipse parameters are 
calculated and subsequent analysis is performed. A problem with this approach is that 
it does not allow cuts on the data to be made in real time based on chi-square test or 
track radius. Thus the output of the ‘low power’ scans tend to large (typically 10,000 to 
100,000 events), because of the number of surface imperfections that axe initially identified 
as track candidates. This has not been considered to be a serious drawback, however, since 
most of the surface imperfections (low energy a-particle tracks, and ‘bubbles’ in the plastic 
surface) are of a spherical or elliptical shape and can not be cut from the data set on the 
basis of chi-square test of the fit alone. Also, the cut on pixel number, while not providing 
as convenient a way to eliminate undesirable features in an image as a cut on minor axis, 
does provide a way to place some cut on object size. The cuts on minor axis can be done 
later on the micro VAX, as the analysis requires. 
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4. Analysis Procedures 


As mentioned earlier, the scan is divided into two parts, a c low power’ scan (typically 
a OX), in which the surface of the plastic is scanned rapidly and rough calculations are 
made of the ellipse parameters, and a ‘high power’ scan (typically 50X), in which the track 
parameters are measured accurately. Given these two magnifications, the difference in the 
accuracy of the determination of ellipse major and minor axes is easily estimated from the 
e( 3ge position resolution at 10X and 50X. At 10X, each pixel corresponds to 2.5 /xm, and 
the number of pixels in a typical track is approximately N = 30. Thus while the position of 
a pixel only specifies a position to ±1.25/im, the position of the edge can be given with an 
uncertainty <7 Edg€ = 1.25 j\fN or 0.23^m. A similar calculation applies for tracks measured 
at SOX. Here the micron/pixel scale factor is approximately 0.5, and the number of pixels 
per track is about N = 400 for most of the samples we have analyzed, so <r Edge = 0.5 fy/N 
or roughly 0.03/xm. While these values for a Edg€ can be realized in practice if successive 
measurements are made while the track is viewed at one position under the camera, the 
actual resolution that we obtain is limited by the effect of image distortion, (±2% over 
the entire field of view). The accuracy in determining the major or minor axis is roughly 
±0.2/xm at SOX. 

The image distortion also must be taken into account when calculating the sheet po- 
sitions (in /im) from the centroids of the objects calculated in pixel coordinates during 
a low power scan. At 10X, at typical low power scan magnification, the field of view is 
1200 X 950/xm~, so that a 2% effect is roughly 20/xm over the entire field of view, a sig- 
nificant fraction of an average track dimension. A linear micron/pixel scale factor is thus 
not sufficient to accurately give the sheet positions of tracks given their centroid x and y 
pixel coordinates. In order to eliminate this effect, we have developed a program to map 
the pixel coordinates to sheet coordinates using a table lookup. This is done by raster 
scanning a single track across the screen and tabulating track centroid pixel coordinates 
for a given stage position. The centroid of the track can be estimated to sub-pixel accu- 
racy (see the arguments above); the sheet coordinates are limited by the accuracy of the 
position measuring apparatus to ±l/im. After this table is accumulated, the results are 
inverted using a bilinear interpolation scheme to give the sheet coordinates for each point 
on a rectangular grid of pixel coordinates. A typical grid dimension is 10 X 10. This grid 
is stored as an array in the analysis program and serves as a look-up table for calculating 
the sheet position of arbitrary pixel values, again using a bilinear interpolation process. 
Using this method we have reduced the effects of camera distortion to the approximately 
.5% level over the viewing area. 

While this look-up table could be used in principle to calculate the sheet coordinates 
for every pixel in a track, this has not been implemented in our system. In general, the 
tracks do not take up a significant area on the image, and, at high power, the tracks 
are well centered so that the micron/pixel scale factor is roughly constant over the area 
of interest, and a scale factor is used to calculate the sheet coordinates of a given pixel. 
Since the scale factors are different for the x and y directions (due to the fact that the 
pixel aspect ratio is 4:3 for the RS-170 television standard), two scale factors must be 
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calculated. This is done in a manner similar to that described above. The average values 
for the x and y pixel/micron conversion factors are calculated for a number of points in 
a region around the center of the image by moving a single track in a raster scan across 
the screen. The position in pixel space of the centroid of the track is compared with the 
stage position, and the average of the micron/pixel scale factor thus deduced is used in 
subsequent calculations. 

The analysis of the plastic begins with a scan at low magnification of both sides of 
a sheet of plastic. This is a raster scan of an entire sheet of plastic (or a region selected 
interactively from the menu), where each ‘frame’ (the image at a particular location in 
the scan) is processed to give sheet coordinates and approximate ellipse parameters for 
all potential tracks identified. Adjacent frames overlap one another by an amount which 
is adjustable from the menu, so that track candidates which lie on the boundary of a 
particular frame will not be excluded. 

The focus is maintained by a Leitz laser autofocus system. For the low power scan, 
the autofocus is not under control of the computer and no analysis time is devoted to 
this task. This gives an improvement in speed over an earlier system which used a ‘focus 
map’ to determine the appropriate focus for each point on the sheet, by interpolating from 
the z-axis position at focus for a grid of (x,y) coordinates on the sheet. Also, this focus 
map had to be constructed each time the sheet was placed on the stage. The autofocus is 
switched on at the beginning of the scan by a control keypad available to the operator. 

The images are first analyzed using a simple gray level threshold to discard images 
which are unlikely to contain tracks. The threshold is selected such that the gray levels of 
the majority of surface defects or non-penetrating tracks, which are shallow and therefore 
not as dark in reflected light as tracks will well developed etch cones, will be below the 
chosen threshold level. The result of this threshold is a binary image (all pixels either 0 
or 1), with the value 1 corresponding to interior points of possible tracks. The number 
of pixels with the value 1 (‘on’ pixels) is then used as a cut to determine whether the 
image requires further processing. This counting is done by histogramming the gray levels 
of the binary image, the contents of the T’ bin in the histogram being the number of 
pixels that correspond to possible track interior points. If the number of ‘on’ pixels after 
thresholding is above a certain value, the image is processed further to segment the edges 
of the track candidates; otherwise, the image is discarded and the next frame in the scan is 
digitized and processed. Because the thresholding and subsequent histogramming can be 
done very fast in hardware (see section 2), this method provides a fast method of skipping 
uninteresting frames in the low magnification scan, if the density of tracks is low enough 
that there are a considerable number of frames without track candidates, and if the surface 
is smooth enough that most of the uninteresting features can be eliminated by a simple 
thresholding process. A cut is also placed on the maximum number of pixels that a frame 
may contain and still be considered for processing. This allows the program to quickly 
skip over areas, such as the mounting holes for the plastic sheets, without subjecting them 
to edge analysis. 

In practice the threshold level and the pixel number cut are determined by viewing 
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the scanned images on the screen and adjusting these parameters interactively so that the 
a gonthm detects most of the events considered interesting by the operator In some ties 

he nt^r ^ * ° aly ^tectLne smdlpo™ 

possibility is t^ Il’a" P “ el “ C “‘ mUSl d ““ 

If the frame passes these pixel number cuts, it is subjected to edge analysis discussed 
previously (see also appendix A). The only real-time cut placed on the ot ects^ 
is a minimum purel number cut, and again this can be set interactively from the menu 
program to select tracks appropriate to the analysis. In this way smaU imTrfSoTs in 
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routine and the seaxch routine run in a time which is proportional to TVlogiV, the entire 
routine runs in essentially TVlogiV time, and a dramatic time saving is achieved over the 
‘brute force’ N~ algorithm. The match is further constrained by placing cuts on major 
and minor axes and azimuth angle, except for tracks with incident angles of approximately 
0 , where the azimuth angle is not specified. Good results have been achieved at 10X by 
requiring that the major and minor axes for the events to be matched agree to within 1/zm, 
and the azimuth angle to within 10 0 , with a search area of 50 X 50^m 2 from the calculated 
point. For these values, more than 75% of the tracks that pass the cuts are ‘good’ events. 

The next phase in the analysis process is the high magnification scan. The files pro- 
duced by the matching algorithm are loaded onto the VICOM via the DMA link and used 
to drive the microscope to the sheet locations of track candidates. The edge segmentation 
operation is essentially the same as that for the low power scan (see section 3). The number 
of tracks to be processed is generally much less than the number processed in a high power 
scan; the track parameters are calculated on the VICOM, and written directly to disk. 
This allows us to perform additional cuts on the measured tracks. As of now the only cut 
that is available is track minor axis, a ‘goodness of fit’ cut has not yet been implemented. 

For the high magnification scan, the autofocus is controlled by the VICOM microcom- 
puter. The autofocus is switched on while the stage is moving to a track coordinate, and 
switched off prior to stopping. This prevents the focus from being lost as the track comes 
into the field of view, and the laser spot is incident on the surface of the track cone. The 
time required for the stage to move from one sheet coordinate to another is easily calcu- 
lated given the velocity profile of the stage. The Compumotor indexers, which send TTL 
control signals to the autofocus via available programmable output bits, can be instructed 
to delay processing of the autofocus commands via am ASCII string sent from the VICOM 
microcomputer. In this way the appropriate delay can be generated before the autofocus 
is disabled upon reaching a track coordinate. 

The time required to complete a high power scan is currently limited by the stage 
movement between successive screen coordinates. The problem of minimizing the travel 
time between a number of points on a surface has been well studied, and an approximate 
solution based on simulated annealing’ is given in Press et al. (1986). Because, as a result 
of our matching method, the high power scan files are sorted according to y coordinate, this 
travel time is certainly not optimized. We plan to implement some sort of solution to this 
problem in the future. Our high magnification scan routine works either automatically, 
or with operator assistance. In the latter case, the movement of the stage to successive 
track locations is initiated by a command sent from the menu task. The operator can also 
choose to analyze a track via another menu call, and also to record the track parameters 
in a file. This method of control has proved useful in debugging our analysis routines. 
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5. System Performance 


We have tested the performance of the track analysis system on CR-39 data obtained 
in accelerator experiments and on a balloon flight. An attempt has been made using the 
accelerator data to determine the charge resolution that can be obtained with the system, 
and to investigate systematic variations in track measurements with sheet position and 
other effects. A full description of these experiments has been written (Dowell 1987), and 
a brief discussion of the results will be given here. The full analysis of the balloon data is 
still in progress, and here a discussion of the results obtained so far will be given. 

The accelerator exposures were done at the Berkeley Bevelac to a a 1.7 A GeV 56 Fe 
beam, using a stack of 6 X6 X 0.072 " CR-39 sheets with an air-gap of 1.13 cm. The stack 
was composed of repeating modules of two sheets in contact, followed by the air gap, for a 
total of 110 sheets. These sheets were etched at 40 ° C in 6.25M NaOH for 480 hours’. The 
resulting Fe primary tracks were approximately 50jma in diameter. A 2" X 2" region of the 
37th sheet in the stack was scanned and approximately 19,000 objects recorded. The scan 
took roughly four days to complete. A cut was placed on the iron peak (R Ft = 24.3/mi), so 
that other peaks could be easily identified. Using the semi-automatic measuring program, 
the fronts of sheets 36 and 38 and the back of sheet 38 were scanned at high magnification, 
using the coordinates from the low power scan of sheet 37. The operator was responsible 
for identifying tracks and eliminating bad events from the data file. These were mainly 
surface imperfections which produced shallow circular pits with a characteristic bright 
spot in the center, and were easily discarded from the data by inspection of the image on 
a video screen. The data from the three surfaces was correlated and and charge changing 
interactions eliminated from the data set. After these cuts the data set consisted of 951 
events. The histograms of these surfaces are shown in figure 5-7. A histogram of average 
radius measurements, figure 8, was used to determine the charge resolution of the CR-39. 
In order to identify the charge peaks in the histograms correctly, a calibration of the Fe 
tracks was made, since the cut on the low magnification data set was severely limited the 
number of Fe and Mn tracks which appeared in the high magnification scans. 

The charge resolution is given by cr/AR, where <r is standard deviation of a single 
charge peak, and A R is the separation between two adjacent peaks in figure 8. The charge 
resolution was calculated for each peak, 12<Z<24, and the mean charge resolution was 
determined to be a z = (0.211 ± 0.0012)/ y/n, where n is the number of scanned surfaces. 
This is about the same resolution reported by other authors (Gerbier, Williams, Price and 
Guoxiao 1987) for CR-39(DOP). Figures 9-11 are plots of Rj vs. i? 2 where Rj and R 2 
represent radii measurements of an individual non-interacting track on two different sheets. 
The individual track measurements differed by less than 1.5/im between two surfaces. 

A study was also made to try to discover any systematic differences in the radii of the 
tracks as a function of sheet position. The sheet was divided into 3" X 3" regions (9 total), 
and the average radius for each charge species in each each region was calculated. The 
weighted distribution of radius variation for each region, normalized to an arbitrary region 
(which we labeled with the number 9) was calculated using: 
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where Rz % is the average radius for charge Z in region i, Rz 9 is the average radius for 
charge Z in region 9, Nz i is the number of radii with charge Z in region i, and N{ tot is the 
total number of tracks in region i. The weighted distributions for each region differed by 
less than 0.5% from region 9. This result implies that variations in track measurements 
with respect to position is negligible. 

The above measurements were done with a system which in which most of the track 
segmentation and fitting routines were written in PASCAL. Since then, in an attempt 
to improve the speed of the analysis, much of the segmentation routine has been written 
in assembly language, and the fitting done on the micro VAX from the raw moment data 
derived on the VICOM. We have also developed the software for matching tracks on ad- 
jacent surfaces of the sheets, and the software for the autofocus to allow us perform high 
power scans automatically. The performance of the improved system has been tested on 
data collected during a balloon flight. 2080 sheets (6" X 6 " X 0.010") of CR-39(DOP) were 
flown from Prince Albert, Saskatchewan, for 5.5 hours at a float altitude of 120,000 feet. 
The plastic was etched for 504 hours in 6.25M NaOH. Subsequent analysis was performed 
on the sheets at low magnification to determine the location and approximate major and 
minor axes of the events. A typical scan provided lists of 60,000 event candidates. Each 
side of the sheet took approximately 24 hours to scan at 10X magnification, this time de- 
pending somewhat on surface quality and number of track candidates found. The matching 
procedure between the front and back surfaces of a single sheet yielded roughly 600 events. 
The cuts placed on the data were: 1.0/im on major and minor axis, and S ° on azimuth 
angle. The search area was set equal to 3 times the major axis of the ellipse. The number 
of events found was not highly dependent on this parameterization, and it was chosen 
because the events at large zenith angle showed larger deviations from calculated positions 
than those with smaller zenith angle. Figure 12 shows a histogram of b (minor axis) for 
one particular sheet. This values for b must be corrected for the deviation in the center of 
the etch cone from the center of the track mouth (a function of zenith angle) to provide 
information on (Z//9). A more convenient way to study this is to plot the ratio of track 
etch rate to bulk etch rate, s = v T /v G which gives a unique value for (Z / 0) independent 
of zenith angle. This is shown in figure 13 for the data of figure 12. The value of s shows 
a sharp peak around the value of 1.5. This corresponds well with the value of s found for 
relativistic iron in a previous study (Tarle et al. 1981). A plot of 8 is shown in figure 14. 
Ideally, the distribution of tracks in theta should have a peak at around 40 ° , while our 
data shows a peak at less than 30 ° . The discrepancy can be attributed to a low detection 
efficiency at large 0, since the tracks at large angle of incidence tend to be less dark than 
those near 8 = 0 and are likely to be missed by the gray level threshold used to pick 
'interesting frames. This is even more likely for tracks with large zenith angle and low 

(Z//J). 



The charge resolution in this experiment is ultimately to be determined by matching 
the events in the plastic with electronic data from the to give positive identification of 
iron events, and measuring the events found at low power with a high magnification scan. 
This work is m progress. Note that a cut on low 0 events is implicitly performed by the 
requirement that the major and minor axes differ by no more than 1.0/im, but the exact 
parameters for this cut have not yet been quantified. 
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Appendix A: Edge Detection Algorithms 

The first stage in the segmentation of an edge concerns the detection of intensity 
changes in an image. The general idea behind this detection is that a sudden intensity 
change in an image will give rise to a peak in the first and second derivative, or a zero 
crossing in the second derivative of an image. Marr and Hildreth (1980) have argued 
that the most satisfactory method for detection of intensity changes is to locate the zero 
crossings in an image convolved with the operator V 2 (? where V 2 is the two dimensioned 
Laplacian operator and G stands for the two dimensional Gaussian distribution 

-f 

G(x,y ) = e 

which has standard deviation <r. V 2 G is a circularly symmetric Mexican hat-shaped opera- 
tor, whose distribution in two dimensions may be expressed in terms of the radial distance 
r from the origin by 



Where I is the raw image data. By the convolution theorem, (V 2 G)*/ = V 2 (G*I), so this 
process can be thought of as the convolution of I with the operator (V 2 G), or equivalently, 
as the operation of V 2 on the smoothed image (G * I). 

There are two basic ideas behind this choice of filter. The first is that the Gaussian 
part of it, G , blurs the image, wiping out all structure on scales much smaller than the 
space constant a of the Gaussian. If the a is chosen small enough, the effect is to eliminate 
much of the noise in the image that is due to small scratches in the plastic, or single 
‘bright' pixels. Marr and Hildreth chose a Gaussian, instead of a pillbox or some other 
smoothing function, because it is a smooth function localized in both the spatial and 
frequency domains, and as such it is unlikely to introduce changes that were not present 
in the original image. 

The second idea concerns the derivative part of the filter, V 2 . In many image processing 
applications, first order derivatives are used to detect edges in some combination like 

)'+(£)’ 

where I(x,y) is the intensity at the point (x,y) in the image, and the partied derivatives 
can be represented by the differences: 

= /(*,y) - /(* - l,y) 

dl 

0 ^ = f{x,y) - f{z,y- 1) 

The above method is computationally very expensive, and some means must be found 
of computing the square and square root operation on an entire image. The following 
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approximation to the above is generally used: 


|v,(*„) W £| + ,£i| 

However this approximation is strongly directionally dependent, and thus some edges with 
Thus'it ° nenlatl °” (fo^bly along the x or y direction) will be enhanced more than others 

edge from background" 11 * threSh °' d Which «*" go ° d res “ its “ «» 

climhndine"d' 0n f “f! “ d P 1 ^* 1 *° ‘ he problem of *»P«Wng the computation and 
op3or7v^Wn‘th^ dePende ““ “ t0 “ Se the spherically symmetric second derivative 
P r ,, ' th a PP r °ach, the edges are then chosen by detecting the zero cross- 

TSkabl' ^ tmg ‘ mage ' This has the advantage over the threshold method in that it is 

men Ah tT/ a ° f COntrasts ' The optimal thr “ h “ld election in an enhance- 
ent / 1 hresholdmg scheme is generaUy sensitive to the edge contrast and if this contrast 

caries widely no one threshold value may be appropriate. On the other hand, the extent 
of a zero crossing does not depend in any way on the contrast of the edge. 

f ®' Kd * m ' thod ° f attracting edges which is based on Marr and Hildreth’s but differs 
m theirs based on the limitations of our hardware and the time allowed for processing 

v' by gxiT '^t Same T ° ! ° PWationS ( V2 °) * 1 = V2 (G ♦ I) but represent G and 
4b d fi d “ n ™lutions, these are shown in figure 4. The smoothing operation figure 

mask wHI “is fh r a ’ ? mateS ° CCUrr!ng OTer these «**>■ The Laplacian 
hw’rh a . i e a PP roJt,malio " to a spherically symmetric function We 

have chosen to detect the maxima of this function, instead of zero crossings since the 

The 'threshold direct' 7 7!,,™?' 7 ° * Wide r “ ge ’ “<* the illumination is uniform, 
ine threshold detection algonthm is described in Appendix B. 

It is interesting to note that a recent theorem of Logan (1977) states that a one octave 
ndpass function can be reconstructed, up to a multiplicative constant, by a knowledge 
of its zero . crossings i alone. While the (V*C) function is certainly not a pure 0^7 lave 
bandpass filter (its bandwidth at half power is 1.25 octaves for the parameters uTed in "he 
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Appendix B. Threshold/Tracking Algorithm 

The t “/ tr - ki ng algorithm that we use to locate edges from the grey level 
maxima m the laplacian image has been briefly described in section 3. Figure 15 shows 

—hi r tit fOT tMS T ^ algorithi prod™ 

,r g h r gat,v r posiiive peaks in the ^ kvd *>■ ^ i°««on of 

traoLr^rihm 1 t . L * P ‘* d “ iJ * orilhm °" i d «d edge. The threshold 

some thrahold T Wh S ^ r ° W , S mem ° ry fOT a negative S™* level below 

some threshold, (. When a pixel which satisfies this criteria has been found, the eight- 

neighbor adjacent pixel, are searched to find a neighbor with the lowest gray levd vflue 

tch S "thf, “V'. ’'I?' " fOUnd ’ U tha *“*« 

search. From the most recent pixel found, 3 neighboring pixels are searched for a grey 
level below the value of t . Figure 16a and 16b show the pixel search list for two given 
directions. Pixel 2 represent the most recent pixel found in the search, and pixel 1 * the 
p xe. found previous to pixel 2. The pixel search list is defined sui that te 3 lei! 

Pixels ltd 2 Th m ° St nearly , fOTmin * a strai S hl Iine ‘he direction defined by the 
pixels 1 and 12 This prevents the search from tracing back over itself, and breaking the 

con mm y o t e track at locations which the gray level minima form ‘blobs’ on the track 

dge. Each successive pixel found in the search defines a new list of pixels to be searched 

”r R gm^ ’ T its locMon “ d diction it defines with the previous 

pixeh Figure 17c shows how this search process proceeds in relation to pixels which make 

flit be , found f Tha r ch continues “" iil a > — 
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bm«T , ^ ? ,S n °‘ * ‘ raCk edga ’ and the list of ‘fo>“ d ’ pixels is deleted 

, n T d 'T n unvranled ed ses found from consideration as possible track 
candidates. One example of this type of edge is that produced by a scratch in the plastic 
which appears as a long straight line of pixels in the processed image. ? 

in tLT^ tha ‘ th A> entire ‘ raCk edge be found may «“ difficult to assure, 

moceX e outh^f h geS " CW ? D0P > « "OAaO, weU defined that the search 
p cedure outlined above gives good results even with tracks produced at large zenith 
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^v urf r: may ; ausa to io - tha ^ “ d »niy a 

H^tvl fmmT l y COn i n b) , aboVe the track is than disc “ded from the analysis, 
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The remainder of this appendix is a code listing for the threshold/ tracking algorithm 
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BUZZ 


1 


i out i, 

NOPAGE 

XDEF BUZZ, PREBUZZ 

SECTION 1 

BUZZ12 incorporates both modifications of BUZZ7 and BUZZ8 . This is 
the algorithm which is currently being used. DC 12/9/88. 

Added HAXPIXEL routine 12/20/88 DC. 

PROCEDURE BUZZ (M INPIXEL, MAXP I XEL : WORD; SPARSE : WORD; VAR NO OBJECTS : WORD 
VAR PIX NO LIST : NO ARRAY; VAR PV : PTRARRAY; - 
VAR RSUH, CSUM : INTARR) ; 

TYPE POINTERl - "PIXEL REC; 

PIXELREC - RECORD 

COORD : ARRAY ( 1 . . 2 ] OF WORD; 

NEXT : POINTERl; 

END; 

PTR ARRAY - ARRAY |1 .. 500 ] OF POINTERl; 

NO ARRAY = ARRAY (1.. 500) OF WORD; 

INTARR - ARRAY ( 1 . . 500] OF INTEGER; 

ROW AND COLUMN VALUES GO FROM 0 TO 511 
CURRENT REGION OF INTEREST : ROWS (5 - 465) 

: COLUMNS (36-455) 

START ADDRESS IS $800000 *( IN ITIAL ROW) *1024+ (INITIAL COLUMN) *2 
BINARY IMAGE IS ASSUMED TO BE IMAGE |0 

PREBUZZ (IN ITROW, INITCOL, FNALROW, FNAL^COL, NUM COL, CCOUNT, RCOUNT : WORD; 

STRT ADD, Zl BEG, Z2 BEG : INTEGER) 


PREBUZZ LOADS THE FOLLOWING MEMORY LOCATIONS: 

NUM COL - FNAL COL-INITCOL 

CCOUNT - FNAL_COL-INIT_COL+2 

RCOUNT - FNAL ROW- INIT ROW 

STRT__ADD - $800000 + INIT ROW* 1 024* IN IT COL*2 

Z1 BEG - $800000* (INITROW-l) *1024*TlNIT COL-1) *2 

Z2BEG - $800000* (INITROW) *1024+ (INITCOL-l) *2 


(PREVIOUS VALUES USED HAVE BEEN: 


* 

PIXELVAL 

EQU 

$01 

INIT ROW 

DC . W 

0 

INIT COL 

DC . W 

0 

FNAL ROW 

DC . W 

0 

FNAL COL 

DC . W 

0 

NUM COL 

DC . W 

0 

CCOUNT 

DC . W 

0 

RCOUNT 

DC . W 

0 

STRT ADD 

DC.L 

0 

Z1 BEG 

DC . L 

0 

Z2 BEG 
* 

DC.L 

0 

STACK1 

DS.L 

200 

STACK2 

* 

DS.L 

200 

HEAPBEG 

DS.L 

50000 

HEAPEND 

EQU 

* 


DC.L 

0 

HSBASE 

DC.L 

0 

ACCUM 

DC . W 

0 

RETPTR 

DC.L 

0 

ARRPTR 

DC.L 

0 


INITROW-5, INIT COL-36 
FNALROW- 465, FNALCOL- 455) 


* RETURN POINTER TO PASCAL 

* POINTER TO ARRAY AREA 



OBJPTR 

DC . L 

0 

NODEPTR 

DC.L 

0 

MINPIXEL 

DC . W 

0 

MAXPIXEL 

DC . W 

0 

EDGFLG 

DC . W 

0 

NUMUPD 

DC . W 

0 

DIRECT 

DC . W 

0 

EDGMARK 

DC . W 

0 

NEWADD 

DC.L 

0 

INITDR 

DC . W 

0 

R SUM 

DC.L 

0 

C SUM 

DC.L 

0 

CLOSED 

DC . W 

0 

DEADFLG 

DC . W 

0 

CHKADD 

DC.L 

0 

JTABLE 

DC.L 

DIR1 


DC.L 

DIR2 


DC.L 

DIR3 


DC.L 

DIR4 


DC.L 

DIRS 


DC.L 

DIR6 


DC.L 

DIR7 


DC.L 

DIRS 


SECTION 

9 

PREBUZZ 

EQU 

* 


MOVE . L 

(A7) ♦ , RETPTR 


MOVE . L 

( A7 ) +, Z2 BEG 


MOVE . L 

(A7) ♦ , Zl BEG 


MOVE . L 

(A7) t, STRT ADD 


MOVE . W 

( A7 ) +, RCOUNT 


MOVE . W 

( A7 ) + , CCOUNT 


MOVE . W 

( A7 ) + , NUM COL 


MOVE . W 

(A7)+,FNAL COL 


MOVE . W 

( A7 ) + , FNAL ROW 


MOVE , W 

( A7 ) +, INIT COL 


MOVE . W 

(A7)+,INIT ROW 


MOVE . L 

RETPTR, A4 


JMP 

(A4) 

BUZZ 

EQU 

* 


* JUMP TABLE FOR SEARCH ROUTINES 


REGISTERS : DO ROW COUNTER 
D1 FINAL COLUMN 

D2 ROW INCREMENT FOR SPARSE SCAN 
D3 COLUMN COUNTER 
D4 MEMORY INCREMENT - <D2*1024) 
D5 LOW ORDER WORD - $0000 


AO POINTER TO CURRENT ROW 
A1 POINTER TO CURRENT PIXEL 


MOVE .L ( A7 ) ♦, RETPTR 

MOVE . L (A7)*,C_SUM 

MOVE . L (A7)+,R SUM 

MOVE L (A7)+,AT * POINTER TO ARRAY OF POINTERS 

MOVE L ( A7) +, ARRPTR * POINTER TO ARRAY OF PIXEL NUMBERS 

MOVE . L ( A7 ) + , OBJPTR * POINTER TO NUMBER OF OBJECTS 

MOVEQ 10, D2 

MOVE . W ( A7 ) ♦, D2 

MOVE . W (A7) +, MAXPIXEL 

MOVE . W (A7) +,MINPIXEL 

MOVEM . L A3/A5-A6, - (A7) * SAVE THESE REGISERS ON THE STACK 

* NOW SET THE BORDER PIXELS OF THE 

* IMAGE TO |4. PRESERVE A4 AND D2 ! 




M,D7 
*1024,06 
Z1_BEG,A0 
CCOUNT, DO 
AO, A1 
D7, (A1 ) + 
DO, LOOP A 

»2,A1 
D6,A1 
A1,A0 
Z2_BEG,A3 
A3, A5 
RCOUNT, DO 

D7, (Al) 

D7, (A5) 
D6,A0 
D6, A3 
A0,A1 
A3,A5 
DO, LOOP B 

CCOUNT, DO 
D7, (A3) -f 
DO,LOOPC 

110,00 
D2,D4 
DO, D4 

ACCOM, A6 
*0, (A6) 
HEAPBEG, A6 
ARRPTR, A5 

STRT ADD, A3 

90, dG 

INIT ROM, DO 

l$FF70,D5 

*0,D3 

10, D1 

10, D6 

FNAL COL, D1 
FNAL ROM, D6 


IIUM COL, D3 
D2 , DO 
06, DO 
DONE 
04, A3 
A3, Al 
A3, A2 

*$80000, A2 
*2, A2 
*0, (Al ) -f 
PROCEED 
-2 (A2) , D5 
CHECK 
D3, LOOP2 
LOOP1 


DONE 


move.l 

MOVE . M 


OB.JPTR, A4 
ACCOM, (A4) 


* D6 


1024 


* D2 - NUMBER OF ROMS SKIPPED 

* D4- (02*1024) 

* CLEAR THE LOCATION FOR NUMBER OF 

* OBJECTS COUNTER 

* A6 IS THE HEAP POINTER 

* AS POINTS TO THE CURRENT 

* POSITION INTHE PIXEL # ARRAY 

* INITIALIZE REGISTERS 


* PIXELVAL - VALUE OF "ON* PIXEL 


* PIXEL NUMBER ARRAY IS OF SIZE MORD 

* D3 IS USED TO COUNT 'INNER' LOOP 

* ADD ROM INCREMENT D2 TO DO 

* ROW + INCREMENT > FINAL ROW? 

* THEN FINISH 

* ADD ADDRESS INCREMENT D4 TO Al 



MOVEM.L (A7) ♦ , A3/A5-A6 

MOVE . L RETPTR,A4 

JMP (A4) 

REGISTERS : DO ROW COUNTER 

D1 COLUMN COUNTER 
D2 COLUMN COUNTER 
D3 JUMP TABLE COUNTER 

AO POINTER TO JUMP TABLE BASE 

Al POINTER TO "CURRENT* PIXEL 

A3 POINTER TO JUMP TABLE ENTRY 
A4 POINTER TO STACK1 

A5 POINTER TO STACK2 

A6 HEAP POINTER 


VO 


SEARCH 


HRESET 


SRCH 


JMPRET 


EQU 

* 



MOVEM.L 

DO-D6/AO-A3, - (A7) 

* 

SAVE EVERYTHING ON THE STACK 

MOVE . L 

C SUM, A3 



CLR.L 

(A3) 

* 

CLEAR ROW AND COLUMN SUM LOCATIONS 

MOVE . L 

R SUM, A3 



CLR.L 

(A3) 



LEA. L 

JTABLE, AO 



MOVE . L 

A7 , HSBASE 

* 

HSBASE-HARDWARE STACK BASE 

SUBQ 

12, Al 

* 

Al IS THE CURRENT PIXEL POINTER 

SUBQ 

• 2 , A2 



SUB.W 

D3,D1 

* 

D 1-COLUMN I OF CURRENT PIXEL 

MOVE . L 

A 1 , NEWADD 



MOVE . L 

Al , CHKADD 



MOVE . W 

D0,-(A7) 



MOVE .H 

D1,-<A7) 



MOVE . L 

A1,-(A7) 



MOVE .If 

18, (Al) 

* 

MARK THIS PIXEL 

ADDQ.L 

•1, (A5) 

* 

INCREMENT NUM OF PIXELS COUNTER 

MOVE . L 

A6, (A4 ) ¥ 

* 

PUSH A NEW NODE POINTER ON 



* 

THE STACK 

MOVE . L 

A6, NODEPTR 

* 

SAVE NODE POINTER IN CASE OF 



* 

DELETE 

MOVE . W 

HO, (A6) + 

* 

PUSH THIS FIRST PIXEL ON THE STACK 

MOVE.lt 

Dl, (A6)+ 



AD DO 

14, A6 



MOVE . L 

AS, -4 (AS) 



BSR 

INIT 

* 

FIND AN INITIAL DIRECTION 

MOVE . W 

•$FFFF , D4 



MOVE . W 

tO, EDGFLG 

* 

RESET EDGE FLAG 

MOVE . H 

DO, DS 

* 

SAVE THE INITIAL ROW/COLUMN IN 

MOVE . K 

Dl , D7 

* 

REGISTERS D6/D7 

MOVE . L 

NEWADD, Al 



MOVE . L 

NEWADD, A2 



ADD . L 

#$80000, A2 



MOVE . W 

DIRECT, D3 



SUBQ . W 

11, D3 



MULU.W 

1 4 , D3 



MOVE . L 

0 (AO, D3 . W) , A3 

* 

JUMP TO SEARCH SUBROUTINE 

JMP 

(A3) 




MOVE . W 

EDGFLG, D 6 

* 

CHECK IF EDGE PIXEL FOUND 

CMP I . W 

11, D6 



BEQ.S 

EDGEH 

* 

BRANCH TO EDGE HANDLER 

MOVE . W 

DIRECT, D6 

* 

CHECK TO MAKE SURE SOME PIXEL WAS 

CMP I .W 

10, D6 

* 

FOUND 

BEQ.S 

DEADEND 

* 

DEAD END HANDLER IF DIRECTION-O 




MOVE . W 

DO, ( A6 ) 4- 

* OTHERWISE WE ARE IN BUSINESS 


MOVE . W 

Dl, (A6) + 

* PUT PIXEL'S ROW, COLUMN COORD. 01 


ADDQ 

14, A6 

* HEAP 


MOVE . L 

A6,-4 ( A6 ) 



MOVE . W 

H, (A3) 

* SET BIT IN CURRENT PIXEL 


MOVE . L 

R SUM, A3 



ADD . L 

0(5 , (A3) 



MOVE . L 

C SUM, A3 



ADD.L 

oT, (A3) 



ADDQ. W 

#1, (A5) 

* INCREMENT THE PIXEL « COUNTER 


ADDQ. If 

t 1 , NUMUPD 


* 

* 

BRA.S 

HRESET 

* ADD THIS CURRENT PIXEL TO STACK 

* 

DEADEND 

CMP A . L 

HSBASE, A 7 



beq 

DISCARD 



MOVE . L 

Al, CHKADD 



MOVE. If 

f 1 , DEADFLG 


* 

* 

BRA.S 

EDGEH2 


* 

EDGEH 

MOVE. If 

t 1 , EDGMARK 



CMP A . L 

HSBASE, A7 

* IS THE STACK EMPTY? 

EDGEH2 

BBQ.S 

SDONE 

* BRANCH - END OF SEARCH 

MOVE .L 

( A7) , NEW ADD 

* POP A CURRENT POINTER OFF STACK 


MOVE. If 

( A7 ) +, Dl 

* 

MOVE .M 

( A7) +, DO 


* 

• 

CMP I . W 

f 4, INITDR 

* FIND A NEW DIRECTION 


BGT.S 

SUBDIR 



ADD I .W 

14, INITDR 



MOVE .If 

INITDR, D3 



MOVE. If 

03, DIRECT 



BRA 

HRESET 


SUBDIR 

SUBI .W 

14, INITDR 



MOVE. If 

INITDR, D3 



MOVE. If 

03, DIRECT 


* 

* 

BRA 

HRESET 


* 

SDONE 

MOVE , L 

10,-4 (A6 ) 

* SEAL OFF THE LINKED LIST 


MOVE. If 

CLOSED, D6 



CMP I . If 

11. DC 



BNE.S 

DISCARD 



MOVE .W 

EDGMARK, D6 

* WAS ONE PIXEL AN EDGE PIXEL? 


CMPI.W 

11, DC 


BEQ. S 

DISCARD 

* BRANCH IF IT WAS 


MOVE . W 

MINPIXEL, D6 

* IS THE PIXEL COUNT OK? 


CMP . W 

(A5),D6 


BHI.S 

DISCARD 

* BRANCH IF IT ISN'T 


MOVE. If 

MAXPIXEL, D6 



CMP.W 

(A5),D6 



BLT.S 

DISCARD 



ADDQ. W 

t 1 , ACCOM 



ADDQ . H 

12, A5 

* INCREMENT AND CLEAR PIXEL NUMBER 


CLR.W 

(A5) 

* COUNTER 


ADD. L 

14, R SUM 



ADD . L 

94, C SUM 



MOVE . L 

R SUM, A3 



CLR.L 

(A3) 




DISCARD 


EDGE 


* 

* 

CLOSE2 

* 

* 

* 

DIR1 


NEW 10 


* 


TEST11 


MOVE . L 

C SUM, A3 

CLR.L 

(A3) 

CLR.W 

CLOSED 

CLR.W 

DEADFLG 

MOVE . L 

HSBASE, A 7 

MOVEM.L 

(A7)+,D0-D6/A0-A3 

BRA 

PROCEED 


MOVE . L 

HODEPTR, A6 

CLR.W 

(A5) 

SUBQ.L 

M,A4 

MOVE . L 

C SUM, A3 

CLR.L 

(A3) 

MOVE . L 

R SUM, A3 

CLR.L 

(A3) 

CLR.W 

EDGMARK 

CLR.W 

DEADFLG 

MOVE . L 

HSBASE, A7 

MOVEM.L 

(A7) D0-D6/A0-A3 

BRA 

PROCEED 


MOVE . W 

I 1 , EDGMARK 

MOVE . W 

| 1 , EDGFLG 

MOVEQ 

10, D3 

BRA 

JMPRET 


MOVE . W 

| 1 , CLOSED 

BRA 

SDONE 


EQU 

* 

MOVE . W 

*0, DIRECT 

MOVEQ 

tO,D5 

CMP.W 

-2 (Al) , D5 

BEQ.S 

NEW10 

MOVEQ 

1 1 , D5 

CMP.W 

-2<A1),D5 

BEQ 

DEADH8 

MOVEQ 

14, D5 

CMP.W 

-2(A1),D5 

BEQ 

EDGE 

MOVEQ 

«8,D5 

CMP.W 

-2 (Al ) , D5 

BEQ 

TST8B 

CMP.W 

-2(A2),D4 

BLT.S 

TEST11 

MOVE . W 

-2 (A2 ) , D4 

LEA 

-2(A1),A3 

MOVE .L 

A3, NEWADD 

MOVE . W 

D6 , DO 

HOVE . W 

D7,D1 

SUBQ.W 

• 1,D1 

MOVE . W 

|8, DIRECT 


MOVEQ 

»0,D5 

CMP.W 

- 102 6 ( Al ) , D5 

BEQ.S 

NEW11 


* REPLACE OLD NODE POINTER 

* CLEAR THE PIXEL • COUNTER 

* RESET NODE POINTER SP 


* IS THIS PIXEL HARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 




MOVEQ 

11, D5 


CMP.W 

- 1026 (Al ) , D5 


BEQ 

DEADH1 


MOVEQ 

14, D5 


CMP.W 

-1026 ( Al ) , D5 


BEQ 

EDGE 


MOVEQ 

10, D5 


CMP.W 

- 1026 (Al ) , D5 


BEQ 

TST81 

NEW11 

CMP.W 

- 1026 (A2 ) , D4 


BLT.S 

TEST12 


HOVE . W 

-1026 (A2) , D4 


LEA 

- 1 026 (Al ) , A3 


MOVE . L 

A3, NEWADD 


MOVE . W 

D6, DO 


MOVE . W 

D7,D1 


SUBQ.W 

II, D1 


SUBQ.W 

11, DO 

* 

MOVE . W 

11, DIRECT 

* 

TEST 12 

MOVEQ 

10, D5 


CMP.W 

-1024 (Al ) , D5 


BEQ. S 

NEW12 


MOVEQ 

11, D5 


CMP.W 

-1024 (Al ) , D5 


BEQ 

DEADH2 


MOVEQ 

14, D5 


CMP.W 

-1024 (Al) ,D5 


BEQ 

EDGE 


MOVEQ 

#8, D5 


CMP.W 

-1024 (Al ) , D5 


BEQ 

TST02 

NEW12 

CMP.W 

-1024 (A2) , D4 


BLE 

JMPRET 


MOVE . W 

-1024 (A2) , D4 


LEA 

- 1024 (Al ) , A3 


MOVE . L 

A3, NEWADD 


MOVE . W 

D6, DO 


MOVE . W 

D7,D1 


SUBQ.W 

11, DO 


MOVE . W 

12, DIRECT 

* 

BRA 

JMPRET 

* 

DIR2 

EQU 

* 


MOVE . W 

10, DIRECT 


MOVE . W 

10, DIRECT 


MOVEQ 

*0,D5 


CMP.W 

-1026 (Al ) , D5 


BEQ. S 

NEW20 


MOVEQ 

» 1 , D5 


CMP.W 

- 1026 (Al ) , D5 


BEQ 

DEADH1 


MOVEQ 

14, D5 


CMP.W 

- 1026 ( Al ) , D5 


BEQ 

EDGE 


MOVEQ 

#8,D5 


CMP.W 

- 1026 ( Al ) , D5 


BEQ 

TST81 

NEW20 

CMP.W 

- 1026 (A2) , D4 


BLT.S 

TEST21 


MOVE . W 

-1026 (A2 ) , D4 


LEA 

- 1026 (Al ) , A3 


MOVE . L 

A3, NEWADD 


MOVE . W 

D6, DO 


* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH- EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 

* RETURN TO SEARCH 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH -EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 




MOVE . W 

D7 , D I 


SUBQ.W 

f I , DO 


SUBQ. W 

1 1 , D1 

* 

MOVE . W 

II, DIRECT 

* 

TEST21 

MOVEQ 

10, D5 


CMP .W 

- 1024 (Al ) , D5 


BEQ . S 

NEW21 


MOVEQ 

11, D5 


CMP . W 

-1024(A1),D5 


BEQ 

DEADH2 


MOVEQ 

14, D5 


CMP .W 

-1024 (Al) ,D5 


BEQ 

EDGE 


MOVEQ 

• 8, D5 


CMP .W 

-1024 (Al) ,D5 


BEQ 

TST02 

NEW2I 

CMP .W 

-1024 (A2) ,D4 


BLT . S 

TEST22 


MOVE . W 

-1024 ( A2 ) , D4 


LEA 

-1024 (Al) , A3 


MOVE . L 

A3, NEWADD 


MOVE . W 

D6, DO 


MOVE . W 

D7, D1 


SUBQ.W 

11, DO 

* 

MOVE . W 

12, DIRECT 

* 

TEST22 

MOVEQ 

1 0 , D5 


CMP .W 

- 1022 ( Al ) , D5 


BEQ.S 

NEW22 


MOVEQ 

1 1 , D5 


CMP .W 

-1022 (Al) ,D5 


BEQ 

DEADH3 


MOVEQ 

#4 , D5 


CMP .W 

-1022 (Al) ,D5 


BEQ 

EDGE 


MOVEQ 

10, D5 


CMP.W 

- 1022 (Al ) , D5 


BEQ 

TST83 

NEW22 

CMP.W 

- 1 022 ( A2 ) , D4 


BLE 

.IMP RET 


MOVE . W 

-1022 ( A2 ) ,D4 


LEA 

- 1022 (Al ) , A3 


MOVE . L 

A3, NEWADD 


MOVE . W 

D6, DO 


MOVE . W 

D7,D1 


SUBQ.W 

• 1,00 


ADDQ.L 

# 1 , Dl 


MOVE . W 

13, DIRECT 

• 

BRA 

JMPRET 

* 

DIR3 

EQU 

* 


MOVE . W 

10, DIRECT 


MOVE . W 

10, DIRECT 


MOVEQ 

10, D5 


CMP .W 

-1024 (Al) ,D5 


BEQ.S 

NEW30 


MOVEQ 

*1,D5 


CMP.W 

- 1 02 4 ( Al ) , D5 


BEQ 

DEADH2 


MOVEQ 

«4,D5 


CMP.W 

-1024 (Al) ,D5 


BEQ 

EDGE 


* ALTER RuW COLUMN COUNTERS 


* IS THIS PIXEL HARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEH CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* RETUPN TO SEARCH 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


♦ BRANCH-EDGE FOUND 



cn 


MOVEQ 

CMP . W 
BEQ 

HEW30 CMP.W 
BLT . S 
MOVE . W 
LEA 

MOVE . L 
MOVE . W 

move.w 

SUBQ.W 

MOVE.W 


TEST31 MOVEQ 
CMP.W 
BEQ. S 
MOVEQ 
CMP.W 
BEQ 
MOVEQ 
CMP.W 
BEQ 
MOVEQ 
CMP.W 
BEQ 

NEW31 CMP.W 
BLT.S 
MOVE.W 
LEA 

MOVE . L 
MOVE.W 
MOVE.W 
SUBQ.W 
ADDQ . W 
MOVE.W 


- 1024 (Al ) , D5 
TST82 

~1024(A2),D4 

TEST31 

-1024 (A2) ,04 
-1024 (A1 ) , A3 
A3, NEWADD 

06, DO 

07, D1 
II, DO 

•2, DIRECT 


#0, 05 

- 1022 ( Al ) , 05 
NEW31 
• 1,05 

- 1 022 ( Al ) , 05 

DEAOH3 

14,05 

-1022 (Al), 05 

EOGE 

18, D5 

- 1022 (Al ) , D5 
TST83 

-1022 (A2 ) , D4 
TEST32 

- 1022 (A2 ) , D4 
-1022{A1),A3 
A3, NEWADD 
D6, DO 
07, D1 
* 1,00 
11,01 
•3, DIRECT 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


# BRANCH -EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 


TEST32 

MOVEQ 

• 0,D5 


CMP.W 

2 (Al) ,D5 


BEQ.S 

NEW32 


MOVEQ 

1 1 , D5 


CMP.W 

2 (Al) ,D5 


BEQ 

DEADH4 


MOVEQ 

• 4,D5 


CMP.W 

2 (Al) ,D5 


BEQ 

EDGE 


MOVEQ 

18, D5 


CMP.W 

2 ( Al ) , D5 


BEQ 

TST84 

NEW32 

CMP.W 

2 ( A2 ) , D4 


BLE 

JMPRET 


MOVE.W 

2 ( A2 ) , D4 


LEA 

2 (Al) , A3 


MOVE.L 

A 3 , NEWADD 


MOVE.W 

06, DO 


MOVE.W 

07, D1 


ADDQ. W 

1 1 , D 1 


MOVE.W 

*4, DIRECT 

* 

BRA 

JMPRET 

* 

DIR4 

EQU 

* 


MOVE.W 

10, DIRECT 


MOVE.W 

10, DIRECT 


MOVEQ 

»0,D5 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH a 1? it HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER ROW COLUMN COUNTERS 

* RETURN TO SEARCH 


* IS THIS PIXEL MARKED IN BINARY 





CMP .W 

(Al ) , 05 

* IMAGE? 



BEQ . S 

NEW 4 0 

* NO? roiMI'RH IT A NEW CANDIDATE 



MOVEQ 

11 , 05 

* MAPPED WITH A 1? IT HAS ALREADY 



CMP.W 

- 1022 ( A 1 ) , 1)5 

* BEEN LOOKED AT 



beq 

DEADH3 




MOVEQ 

14, D5 




CMP . W 

- 1 022 ( A1 ) , 05 




beq 

EDGE 

* BRANCH-EDGE FOUND 



MOVEQ 

18, D5 




CMP .W 

-1022 (Al) ,D5 




BEQ 

TST83 

| 


NEW40 

CMP . W 

“1022 (A2) ,D4 

* NEW MINIMUM VALUE 



BLT.S 

TEST41 

* NEW ADDRESS IN BINARY IMAGE 



MOVE . W 

-1022 (A2) , D4 

* NEW ADDRESS IN LAP LAC I AN IMAGE 



LEA 

-1022 (Al) , A3 




MOVE. L 

A3, NEWADD 




MOVE . W 

06, DO 




MOVE . W 

07, D1 




SUBQ.W 

11, DO 




ADDQ . W 

11, D1 



* 

MOVE . W 

1 3, DIRECT 



* 

TEST 4 1 

MOVEQ 

»0, D5 

* IS THIS PIXEL MARKED IN BINARY 



CMP .W 

2 (Al) ,D5 

* IMAGE 



BEQ . S 

NEW41 

* NO? CONSIDER IT A NEW CANDIDATE 



MOVEQ 

#1, D5 

• MARKED WITH A 1? IT HAS ALREADY 



CMP .W 

2 (Al) ,D5 

• BEEN LOOKED AT 



beq 

DEADH4 




MOVEQ 

14, D5 




CMP .If 

2 ( Al ) , D5 




BEQ 

EDGE 

* BRANCH-EDGE FOUND 



MOVEQ 

»e, D5 


55 


CMP .W 

2 ( Al ) , D5 




BEQ 

TST84 



NEW41 

CMP.W 

2 (A2) ,D4 

* NEW MINIMUM VALUE 



BLT.S 

TEST42 

* NEW ADDRESS IN BINARY IMAGE 



MOVE .W 

2 ( A2 ) ,P4 

* NEW ADDRESS IN LAP LAC I AN IMAGE 



LEA 

2 (Al) , A3 




MOVE . L 

A3, NEWADD 




MOVE . W 

I>6, DO 




MOVE. If 

07, D1 




ADDQ.W 

11, D1 

* ALTER ROW COLUMN COUNTERS 


* 

MOVE . W 

» 4, DIRECT 



* 

TEST 4 2 

MOVEQ 

»0,D5 

* IS THIS PIXEL MARKED IN BINARY 



CMP.W 

1026 (Al ) ,D5 

* IMAGE 



BEQ.S 

NEW42 

* NO? CONSIDER IT A NEW CANDIDATE 



MOVEQ 

»l,D5 

* MARKED WITH A 1? IT HAS ALREADY 



CMP.W 

1 026 ( A 1 ) ,D5 

* BEEN LOOKED AT 



BEQ 

DEADH5 




MOVEQ 

14, D5 




CMP.W 

1026 (Al ) , D5 




BEQ 

EDGE 

* BRANCH -EDGE FOUND 



MOVEQ 

19, D5 




CMP.W 

1026 (Al ) , 05 




BEQ 

TST85 



NEW42 

CMP.W 

1026 (A2 ) , D4 

* NEW MINIMUM VALUE 



BLE 

iMPRET 

* NEW ADDRESS IN BINARY IMAGE 



MOVE . W 

1 02 6 (A2 ) , D4 

* NEW ADDRESS IN LAPLACI AN IMAGE 



LEA 

1 026 ( Al ) ,A3 




MOVE . L 

A3, NEWADD 




MOVE . W 

06, DO 




MOVE . W 

07 , Dl 




ADDQ.W 

11, DO 




DIR5 


NEW50 


TEST51 


cn 


NEW51 


TEST52 


adpu . w 

» 1 , 1)1 

* 

ALTER Pa 'W (7*1.0(411 COUNTERS 

MOVE . W 

|5, DIRECT 



BRA 

JMPRET 

« 

RETOi.-n TO SEARCH 

EQU 

* 



MOVE . W 

10, DIRECT 



MOVE . W 

10, DIRECT 



MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP.W 

2(A1),D5 

* 

IMAGE 

BEQ.S 

NEW50 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

»1,D5 

* 

MARKED WITH A 1 ? IT HAS ALREADY 

CMP .W 

2(A1) ,P5 

* 

BEEN LOOKED AT 

BEQ 

DEADH4 



MOVEQ 

• 4,D5 



CMP.W 

2 (Al) ,D5 



BEQ 

EDGE 

* 

BRANCH-EDGE FOUND 

MOVEQ 

• 8,D5 



CMP. If 

2 (Al) ,D5 



BEQ 

TST04 



CMP . W 

2(A2),D4 

* 

NEW MINIMUM VALUE 

BLT.S 

TEST51 

* 

NEW ADDRESS IN BINARY IMAGE 

MOVE . W 

2 (A2 ) , D4 

* 

NEW ADDRESS IN LAP LAC I AN IMAGE 

LEA 

2 (Al) , A3 


MOVE . L 

A3, HEHADD 



MOVE. W 

D6, DO 



MOVE . W 

07, D1 



ADDQ.W 

• l.DI 



MOVE . W 

» 4, DIRECT 




MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP .W 

1026 (Al) ,D5 

* 

IMAGE 


BEQ.S 

NEWS 1 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

II, D5 

* 

MARKED WITH 

A 1? IT HAS ALREADY 

CMP.W 

1 02 6 ( Al ) , D5 

* 

BEEN LOOKED 

AT 

BEQ 

DEADH5 




MOVEQ 

#4,D5 




CMP .W 

1 02 6 (Al ) ,DS 




BEQ 

EDGE 

* 

BRANCH-EDGE 

FOUND 

MOVEQ 

#8,05 




CMP .W 

1026 (Al ) , D5 




BEQ 

TST0 5 




CMP.W 

1 026 ( A2 ) , D4 

* 

NEW MINIMUM 

VALUE 

BLT.S 

TEST52 

* 

NEW ADDRESS 

IN BINARY IMAGE 

HOVE . W 

1026 (A2) , D4 

4 

NEW ADDRESS 

IN LAP LAC IAN IMAGE 

LEA 

1026(A)) , A3 




MOVE . L 

A3, NEWAPD 




MOVE . W 

D6, DO 




MOVE . W 

D7, PI 




ADDQ . W 

# 1 , DO 

* 

ALTER ROW COLUMN COUNTERS 

ADDQ.W 

• 1,01 




MOVE . W 

15, DIRECT 





MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP.W 

1024(A1),D5 

* 

IMAGE 

BEQ.S 

NEW52 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

11, D5 

* 

MAPPED WITH A 1? IT HAS ALREADY 

CMP.W 

1024 (Al) ,D5 

* 

BEEN LOOKED AT 

BEQ 

r EADH6 



MOVEQ 

• 4 , D5 



CMP.W 

1024 (Al ) ,D5 



BEQ 

EDGE 

* 

BRANCH-EDGE FOUND 

MOVEQ 

18, D5 



CMP.W 

1024 (Al) , D5 



BEQ 

TST86 





NEW 5 2 


cn 


0IR6 


NEW60 


TEST61 


NEW61 


TEST62 


CMP . W 

1H24 ( A2 ) , D 4 

* 

NEW MINIMUM VALUE 

RLE 

IMP PET 

* 

HEW ADDRESS IN B I NAP Y IMAGE 

move . w 

1024 ( A2 ) , D4 

* 

NEW ADDRESS IN LAP LAC I AN IMAGE 

LEA 

1024 ( A 1 ) , A3 



move.l 

A3, NEWADD 



MOVE. W 

D6, DO 



MOVE . W 

D7,D1 



ADD. W 

11, DO 

* 

ALTER ROW COLUMN COUNTERS 

MOVE . W 

16, DIRECT 



BRA 

JMPRET 

* 

RETURN TO SEARCH 

EQU 

* 



MOVE . W 

•0, DIRECT 



MOVE . W 

10, DIRECT 



MOVEQ 

10, 05 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP .W 

102 6 ( Al ) , D5 

* 

IMAGE 

BEQ. S 

NEW60 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

1 1 , D5 

* 

MARKED WITH A 1? IT HAS ALREADY 

CMP .W 

1026 (Al) , D5 

* 

BEEN LOOKED AT 

BEQ 

DEADH5 



MOVEQ 

#4 , D5 



CMP.W 

1026 (Al ) ,D5 



BEQ 

EDGE 

* 

BRANCIUEDGE FOUND 

MOVEQ 

18, D5 



CMP.W 

1026 (Al ) , D5 



BEQ 

TST85 



CMP.W 

1 026 ( A2 ) , D4 

« 

NEW MINIMUM VALUE 

BLT.S 

TEST61 

* 

NEW ADDRESS IN BINARY IMAGE 

MOVE . W 

1026 (A2 ) ,D4 

* 

NEW ADDRESS IN LAP LAC I AN IMAGE 

I.EA 

1 026 ( Al ) , A3 



MOVE.L 

A3, NEWADD 



MOVE . W 

D6, DO 



MOVE . W 

D?,D1 



ADDQ.W 

11, DO 



ADDQ . W 

» 1 , Dl 



MOVE . W 

15, DIRECT 




MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP.W 

1024 (Al ) , D5 

* 

IMAGE 


BEQ.S 

NEW61 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

11, D5 

* 

MARKED WITH 

A 1? IT HAS ALREADY 

CMP.W 

1024 (Al ) , D5 

* 

BEEN LOOKED 

AT 

BEQ 

DEADH6 




MOVEQ 

14, D5 




CMP.W 

1024 (Al) ,D5 




BEQ 

EDGE 

* 

BRANCH EDGE 

FOUND 

MOVEQ 

*8, D5 




CMP.W 

1024 (Al ) , D5 




BEQ 

TST86 




CMP.W 

1024 (A2) ,D4 

* 

NEW MINIMUM 

VALUE 

BLT . S 

TEST62 

* 

NEW ADDRESS 

IN BINARY IMAGE 

MOVE . W 

1024 (A2) , D4 

« 

HEW ADDRESS 

IN LAP LAC I AN IMAGE 

LEA 

1024 (Al) , A3 




MOVE . I. 

A3, NEWADD 




MOVE . W 

D6, DO 




MOVE . W 

P7, Dl 




ADD.W 

11, DO 

* 

ALTER ROW COLUMN COUNTERS 

MOVE . W 

16, DIRECT 





MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN BINARY 

CMP.W 

1022 (Al) , D5 

* 

IMAGE 

BEQ.S 

NEW62 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

f 1,D5 

♦ 

MARKED WITH A 1 ? IT HAS ALREADY 

CMP.W 

1022 (Al ) ,D5 

* 

BEEN LOOKED AT 



BP.ANv H- EDGE FOUND 


**KU l'KAI'117 

MOVF.Q 14, [‘5 

CMK.W I022<AI),D5 

BEQ FOGE 

MOVEQ #8,D5 

CMP.W 1 022 (A 1 ) , D5 

BEQ TST87 

MEW62 CMP.W l 022 (A2 ) , D4 

BLE JMPRET 

MOVE . W 1022<A2),D« 

LEA 1022 (A1 ) , A3 

MOVE.L A3, NEWADD 

MOVE . W D6, DO 

MOVE.W D7,D1 

ADDQ.W fl, DO 

SUBQ.W | 1 , Dl 

MOVE.W #7, DIRECT 

BRA JMPRET 

« 

D l R 7 EQU * 

MOVE.W 10, DIRECT 

MOVE.W |0, DIRECT 

MOVEQ 10, D5 

CMP.W 1024 (Ai ) , D5 

BEQ.S NEM70 

MOVEQ #1,05 

CMP.W 1024 (Al ) , D5 

BEQ DEADH6 

MOVEQ #4, D5 

CMP.W 1024 (Al ) , D5 

BEQ EDGE 

MOVEQ #0,D5 

CMP.W 1024 (AI ) , D5 

cn BEQ TST06 

CO MEW /O CMP.W 1024 (A2),D4 

BLT.S TEST 7 1 

MOVE.W 1024 (A2) ,D4 

LEA 1 024 ( Al ) , A3 

MOVE.L, A3, NEWADD 

MOVE.W D6 , DO 

MOVE.W 07, Dl 

ADDQ.W 1 1 , DO 

MOVE.W 16, DIRECT 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAPLACIAN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* RETURN TO SEARCH 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1 7 IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAPLACIAN IMAGE 


TEST? 1 MOVEQ «0,D5 

CMP.W 1022 (Al) , 1*5 

BEQ.S NEW71 

MOVEQ #1,05 

CMP.W 1 022 ( Al ) , D5 

BEQ DEADH7 

MOVEQ # 4 , D5 

CMP.W 1 022 (Al ) , D5 

beq edge 

MOVEQ # 8 , D5 

CMP.W 1022 (Al ) , D5 

beq TST07 

NEW? 1 CMP.W 1022(A2),D4 

BLT.S TEST72 

MOVE.W 1022(A2) , D4 

LEA 1 022 (Al ) , A3 

MOVE . L A 3 , NEWADD 

MOVE.W 1*6, DO 

MOVE.W D7,I)t 

ADDQ.W | 1 , DO 

SUBQ.W # 1 , Dl 


* IS THIS PIXEL MARKED IN BINARY 

* IMAGE 

* NO? CONSIDER IT A NEW CANDIDATE 

* MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAPLACIAN IMAGE 


* ALTER ROW COLUMN COUNTERS 



TEST72 


cn 


NEW7 2 


DIR8 


NEWRO 


TEST81 


MOVE . W 

« ; , D I RE<T 



MOVEQ 

10, D5 

♦ 

IS THIS PIXEL MARKED IN BINARY 

CMP . W 

-2<A1),D5 

* 

IMAGE 

BEQ.S 

NEW7 2 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

11, D5 

* 

MARKED WITH A 1? IT HAS ALREADY 

CMP .W 

-2 (Al) , D5 

* 

BEEN LOOKED AT 

BEQ 

DEADH0 



MOVEQ 

f 4 , D5 



CMP .W 

-2 (Al ) , 05 



BEQ 

EDGE 

* 

BRANCH -EDGE FOUND 

MOVEQ 

• 8 , D5 



CMP.W 

-2 (Al) ,D5 



BEQ 

TST00 



CMP.W 

-2 (A2 ) , D4 

* 

NEW MINIMUM VALUE 

BLE 

JMPRET 


NEW ADDRESS IN BINARY IMAGE 

MOVE . W 

-2 (A2 ) , D4 

* 

NEW ADDRESS IN LAP LAC I AN IMAGE 

LEA 

-2 (Al) , A3 



MOVE . I. 

A3, NEWADD 



MOVE . W 

D6, DO 



MOVE . W 

D7, D1 



SUBQ.W 

#1,D1 

* 

ALTER ROW COLUMN COUNTERS 

MOVE . W 

V 8 , DIRECT 



BRA 

JMPRET 

« 

RETURN TO SEARCH 

F.QU 

* 



MOVE . W 

10, DIRECT 



MOVE . W 

10, DIRECT 



MOVEQ 

10, D5 

• 

IS THIS PIXEL MARKED IN BINARY 

CMP.W 

1022 (Al) ,D5 

* 

IMAGE 

BEQ.S 

NEW0O 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

# 1 , D5 

* 

MARKED WITH A 1? IT HAS ALREADY 

CMP.W 

1022 (Al) ,D5 

* 

BEEN LOOKED AT 

BEQ 

DEADH7 



MOVEQ 

14, D5 



CMP.W 

1022 (Al) , D5 



BEQ 

EDGE 

* 

BRANCH=EDGE FOUND 

MOVEQ 

10, D5 



CMP.W 

1 022 ( Al ) , D5 



BEQ 

TST8 7 



CMP.W 

1022 (A2) , D4 

* 

NEW MINIMUM VALUE 

BI>T . S 

TEST81 

* 

NEW ADDRESS IN BINARY IMAGE 

MOVE . W 

1 02 2 ( A2 ) , D4 

* 

NEW ADDRESS IN LAP LAC I AN IMAGE 

LEA 

1022 (Al) , A3 



MOVE . L 

Al, NEWADD 



MOVE . W 

D6, DO 



MOVE . W 

07, D1 



ADDQ.W 

1 1 , DO 



SUBQ.W 

»1,D1 



MOVE . W 

17, DIRECT 




MOVEQ 

10, D5 

* 

IS THIS PIXEL MARKED IN 

BINARY 

CMP.W 

-2 (Al) ,D5 

* 

IMAGE 


BEQ.S 

NEWS 1 

* 

NO? CONSIDER IT A NEW CANDIDATE 

MOVEQ 

1 1 , D5 

* 

MARKED WITH A 1? IT HAS 

ALREADY 

CMP .W 

? (Al ) , 1)5 

* 

BEEN LOOKED AT 


BEQ 

DEADH8 




MOVEQ 

»4,D5 




CMF . W 

-2 (Al) ,05 




BEQ 

EDGE 

* 

BP AIR 71 - EDGE FOUND 


MOVEQ 

10, D5 




CMP.W 

-2 (Al) , D5 




BEQ 

T5T80 




CMP.W 

-2 (A2) , D4 


NEW MINIMUM VALUE 



riEwni 



BLT.S TESTS 2 

MOVE . W - 2 (A? ) , [)4 

LEA -2<A1),A3 

MOVE.L A3 , NEWADD 

MOVE . W D6 , DO 

MOVE . W D7 , D 1 

SUBQ.W II, DI 

MOVE . W 18, DIRECT 


TEST82 MOVEQ 10, D5 

CMP.W - 1 02 6 ( A 1 ) , D5 

BEQ.S NEW82 

MOVEQ 9 1 , D5 

CMP.M - 1 026 (A1 ) , D5 

beq DEADH 1 

MOVEQ 14, D5 

CMP.W -1026 (Al ) , D5 

beq edge 

MOVEQ |0, D5 

CMP.W -1026 (Al ) , D5 

BEQ TST81 

NEH02 CMP.W - 1026 (A2 ) , D4 
BLE TMPRET 

MOVE . W - 1 026 (A2 ) , D4 

LEA -1026 (Al ) , A3 

MOVE.L A3, NEWADD 

MOVE . W 06, DO 

MOVE . W 07, D1 

SUBQ.W | I , DO 

SUBQ.W' fl, D1 

MOVE . W I 1 , DIRECT 

O BRA 3MPRET 


IN IT 


INIT1 


1NIT2 


1NIT3 


1 N IT 4 


MOVE.L |0, D3 

MOVE . W I $FFC0, 04 

MOVE.L NEWADD, Al 

MOVE.L A 1 , A2 

ADD.L 9 $00000, A2 

CHPI.W 90, - 1026 (Al ) 

BNE.S INIT2 

CMP.W - 1026 (A2 ) , P4 

BI.E.S IN1T2 

MOVE . W - 1026 (A? ) , D4 

MOVE . If I 1 , DIRECT 

CMP I . W 10, - 102 4 ( A 1 ) 

BNE.S INIT3 

CMP.W -1024 ( A2 ) , D4 

BLE.S INIT3 

MOVE.W -1024 (A2) ,D4 

MOVE.W 12, DIPECT 

CMPI.W 10, - 1022 (Al ) 

BNE.S INIT4 

CMP.W 1022(A2) ,D4 

BLE.S 1NIT4 

MOVE.W -1022(A2),D4 

MOVE.W 93, DIRECT 

CMPI.W 1 0 , 2 ( A 1 ) 

BNE.S INITS 

CMP.W 2(A2),D4 

BLE.S IHIT5 


NEW ADDRESS IN BINARY IMAGE 
* NEW ADDRESS IN LAP LAC I AN IMAGE 


* ALTER Row COLUMN COUNTERS 


* 13 THIS PIXEL MARKED IN BINARY 

* IMAGE 

I I !? 7 COMSIDER IT A NEW CANDIDATE 
MARKED WITH A 1? IT HAS ALREADY 

* BEEN LOOKED AT 


* BRANCH-EDGE FOUND 


* NEW MINIMUM VALUE 

* NEW ADDRESS IN BINARY IMAGE 

* NEW ADDRESS IN LAPLACIAN IMAGE 


* ALTER ROW COLUMN COUNTERS 


* RETURN TO SEARCH 



INI 1 b 


INIT6 


INIT7 


INIT0 


CHOOSE 


CHECK 


DEADH1 


MOVE . W 

<' (A2) ,D4 

MOVE . W 

14, DIRECT 

CMP I . W 

|0, 1026 ( A 1 ) 

BNE . S 

IN1T6 

CMP .W 

1026 ( A2 ) , D4 

BLE . S 

INIT6 

MOVE . W 

1026 ( A2 ) , D4 

MOVE . W 

•5, DIRECT 

CMP I .H 

10, 1024 (Al) 

BNE . S 

INIT7 

CMP .W 

1024 ( A2 ) , D4 

BLE.S 

IN1T7 

MOVE . W 

1 02 4 ( A2 ) , D4 

MOVE . W 

16, DIRECT 

CMP I . M 

10, 1022 (Al) 

BNE . S 

INIT8 

CMP.W 

1 022 (A2 ) , D4 

BLE.S 

INIT0 

MOVE . W 

1022 (A2 ) , D4 

MOVE . M 

17, DIRECT 

CMP I . W 

• 0, -2 (Al) 

BNE.S 

CHOOSE 

CMP.W 

-2(A2),D4 

BLE.S 

CHOOSE 

MOVE . W 

-2(A2),D4 

MOVE. If 

18, DIRECT 

CMP I . W 

10, DIRECT 

BEQ 

SDONE 

MOVE . W 

DIRECT, D3 

MOVE . W 

D3, INITDR 

RTS 



equ 

* 

CMP l . W 

»0, -8 ( Al ) 

BNE 

PROCEED 

CMP I . W 

|0,-6(A1) 

BNE 

PROCEED 

CMP I .W 

*0, - 4 ( Al ) 

BNE 

PROCEED 

CMP I . W 

#0, (Al) 

BNE 

PROCEED 

CMP I .W 

•0, 2(A1) 

BNE 

PROCEED 

CMP I . W 

10,4 (Al) 

BNE 

PROCEED 

CMP I . W 

10, -2050 (Al ) 

BNE 

PROCEED 

CMP I . W 

10,-1026 (Al ) 

BNE 

PROCEED 

CMP I .W 

«0, 1022 (Al) 

BNE 

PROCEED 

CMP I .W 

1 0 , 204 6 ( A 1 ) 

BNE 

PROCEED 

BRA 

SEARCH 


BSR 

CHK8 

CMP.W 

1 1 , DEADFLG 

BNE 

DISCARD 

LEA 

- J026(A1) , A2 

CMP.L 

CHKADD, A2 

BNE 

DISCARD 

MOVE . W 

1 1 , CLOSED 

BRA 

SDONE 



DEA0H2 BSR 

CMP.W 

BNE 

LEA 

CMP.L 

BNE 

MOVE . W 
BRA 


CHK8 

I 1 , DEADFLG 
DISCARD 
-1024 (A1 ) , A2 
CHKADD, A2 
DISCARD 
I 1 , CLOSED 
SDONE 


CD 

DO 


DEADH3 BSR 

CMP .M 

BNE 

LEA 

CMP.L 

BNE 

MOVE . W 
BRA 


DEADH4 BSR 

CMP .W 

BNE 

LEA 

CMP.L 

BNE 

MOVE . W 
BRA 


CHK8 

#1, DEADFLG 
DISCARD 
- 1022 ( A 1 ) , A2 
CHKADD, A2 
DISCARD 
# 1 , CLOSED 
SDONE 


CHK8 

< I, DEADFLG 
DISCARD 
2 < A1 ) , A2 
CHKADD, A2 
DISCARO 
1 1 , CLOSED 
SDONE 


DEADH5 BSR 

CMP .W 
BNE 
LEA 
CMP . L 
BNE 

MOVE. W 
BRA 


* 


< HK8 

f 1 , DEADFLG 
DISCARD 
1026 (A1 ) , A2 
CHKADD, A2 
DISCARD 
• 1 , CLOSED 
SDONE 


DEADH6 BSR.S 
CMP.W 
BNE 
LEA 
CMP.L 
BNE 

MOVE . W 
BRA 


DEADII7 BSR.S 

CMP .W 
BNE 
LEA 
CMP . L 
BNE 

move . w 

BRA 


OHK8 

1 1 , DEADFLG 
DISCARD 
1 02 4 ( A1 ) , A2 
r HKADD, A2 
DISCARD 
# 1 , CLOSED 
SDONE 


< HK8 

• 1 , DEADFLG 

DISCARD 

1022(A1),A2 

CHKADD, A2 

DISCARD 

# 1 , CLOSED 

SDONE 



* 

DEADH9 

BSR . S 

CHK8 


CMP . W 

ff 1 , DEADFLG 


BNE 

DISCARD 


LEA 

-2 (Al) , A2 


CMP . I, 

CHKADD, A2 


BHE 

DISCARD 


MOVE . W 

f 1 , CLOSED 

ft 

ft 

BRA 

SDONE 

ft 

CHK8 

CMP I .W 

1 8 , - 1 02 6 (Al ) 


BNE.S 

CHK8A 


LEA 

- 1026 (Al ) , A2 


CMP.L 

CHKADD, A2 


BEQ 

RTS 

CHK0R 

CHK0A 

CMP I . W 

18, -1024 (Al) 


BNE.S 

CHK8B 


LEA 

-1024 (Al ) , A2 


CMP .L 

CHKADD, A2 


BEQ 

RTS 

CHK8R 

CHK0B 

CMP I .W 

10, - 1022 (Al ) 


BHE . S 

CHK8C 


LEA 

-1022(A1),A2 


CMP.L 

CHKADD, A2 


BEQ 

RTS 

CHK8R 

CHK8C 

CMP I . W 

1 8 , 2 (Al) 


BNE.S 

CHK0D 


LEA 

2 (Al) , A2 


CMP.L 

CHKADD, A2 


BEQ 

RTS 

CHK8R 

CHK0D 

CMP I . W 

10, 1026 (Al) 


BNE.S 

CHK8E 


LEA 

1026 (Al ) , A2 


CMP . L 

CHKADD, A2 


BEQ 

RTS 

CHK8R 

CHK8E 

CMP I ,W 

18, 1024 (Al) 


BNE . S 

CHK8F 


LEA 

1024 (Al ) , A2 


CMP . L 

UK ADD, A2 


BEQ 

RTS 

CHK0R 

CHK8E 

CMP I . W 

18, 1022 (Al) 


BNE.S 

CHK8G 


LEA 

1 022 ( Al ) , A2 


CMP . L 

CHKADD, A2 


BEQ 

RTS 

CHK8R 

CHK8G 

CMP I . W 

18,-2 (Al) 


BNE.S 

CHK0H 


LEA 

2(A1),A2 


CMP . L 

CHKADD f A2 


BEQ 

CHK8R 

CHK8II 

RTS 


CHK8R 

ADDQ.L 

14, A7 


CMP I . W 

1 1 , DEADELG 


BEQ 

DISCARD 


BRA 

CLOSE2 


RESET SYSTEM STACK 



TST81 

cmp I . W 

f 1 , DEADFLG 


BEQ 

DISCARD 


LEA 

- 1026 (Al ) ,1 


CM PA . L 

CHKADD, A3 


BHE 

DISCARD 


BRA 

CLOSE2 

TST82 

CMP I .W 

1 1 , DEADFLG 


BEQ 

DISCARD 


LEA 

-1024 ( A 1 ) t I 


CMP A . L 

CHKADD, A3 


BNE 

DISCARD 


BRA 

CLOSE2 

TST03 

CMP I . W 

f 1, DEADFLG 


BEQ 

DISCARD 


LEA 

-1022 (Al) , P 


CMPA.L 

CHKADD, A3 


BHE 

DISCARD 


BRA 

CLOSE2 

TST84 

CMPI.W 

1 1 , DEADFLG 


BEQ 

DISCARD 


LEA 

2 (A1 ) , A3 


CMPA.L 

CHKADD, A3 


BHE 

DISCARD 


BRA 

CLOSE2 

TSTB5 

CMP I . W 

1 1 , DEADFLG 


BEQ 

DISCARD 


LEA 

1026 (Al ) , A3 


CMPA.L 

CHKADD, A3 


BNE 

DISCARD 


BRA 

CLOSE2 

TST06 

CMPI.W 

f 1, DEADFLG 


BEQ 

DISCARD 


LEA 

1024 (Al ) , A3 


CMPA.L 

CHKADD, A3 


BHE 

DISCARD 


BRA 

CLOSE2 

TST8 7 

CMP T .W 

f 1 , DEADFLG 


BEQ 

DISCARD 


LEA 

1 022 ( Al ) , A3 


CMPA.L 

CHKADD, A3 


BHE 

DISCARD 


BRA 

CLOSE2 

TST88 

CMP I . W 

11 , DEADFLG 


BEQ 

DISCARD 


LEA 

-2 (Al) , A3 


CMPA.L 

CHKADD, A3 


BNE 

DISCARD 


BRA 

END 

CLOSE 2 



Appendix C: Ellipse Fitting 

The algorithm we use for fitting an ellipse to our list of data points is discussed in a 
paper by Agin (1980). The most general equation for an ellipse is given by 

Ax 2 4- Bxy + Cy 2 + Dx + Ey + F = 0. 

For a given pixel coordinate (x,,y,) we let the pointwise error be given by 

y,) = Ax 2 4- Bxiyi 4- Cy 2 4- Dx{ 4- Eyi 4- F. 

Since this quantity is not positive definite, it must be squared to give the aggregate error, 

e = ^2 £ 2 = + Cyl + Dx i + E y* + E) 2 

where the sum runs over all points in the ellipse. The idea here is to minimize this aggregate 
error by chosing appropriate values of the coefficients (.4, B , . . .). Tins is done by obtaining 
the partial derivatives with respect to the coefficients and setting the resulting equations 
equal to zero. The result is five equations in the five unknown coefficients: 

AM 04 4- BMn + CM 22 + DM 30 4- EM 21 4- M 20 = 0 

AA /31 4- BM 22 4- CM \3 4" DM 21 4" EM \2 4" Afn = 0 

AM22 + BMu + CM 04 4- DM\2 4- EMqj + Mq 2 — 0 

AM30 4 " B I > ] 4 " CMu 4 " D AI20 + EM\\ 4 " A/10 — 9 

. 4^/21 + BM \2 4" CA /03 4- DM\\ 4" EM 02 4" il /01 = 0 

where M 04 = ^31 = £xfy,-, etc. The general solution to this problem has been 

solved on the symbolic manipulation program REDUCE, and the solution is given at the 
end of this appendix. The solution requires 15 moments to be summed and for this reason 
we refer to it as the “15 moment” solution. 

In view of the complexity of the solution above, we have chosen to sum the ellipse 
moments in a coordinate system which has its origin at the centroid of the track of interest. 
The result is that the coefficients D and E in the equation for the ellipse above can be set 
equal to zero, and the equation for an ellipse becomes 


.4r 2 4- Bxy 4- Cy 2 — 1=0 


where we have set F = —1 so that .4 = ( 1/a*) and B = (1/6 2 ) for an ellipse with semi-major 
axis parallel to the x-axis, this is the standard form for and ellipse with the semi-major 
axis and semi-minor axes given by a and 6, respectively. 
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The aggregate error then becomes 


“ - ~ + Bx iy% + Cy }) 2 

and the minimization results in three equations in three unknowns 

AJl/40 + BM$\ + CM22 ^ A/20 


AM 2 j + BM22 + CM U = M u 
AM22 + 5A/]3 + CM04 = A/02 

melt C “ eas “ y . be i s ° lved , ^‘ically for the unknown coefficients. The number of mo- 
ments required for the solution to this problem is 10, the 8 moments in the equations 

centroid' This' "f "j “° men ‘ s (A/|I) M<l M 0i), used in calculating the ellipse 

centroid -This is referred to as the “10 moment solution”. The solution for A B and C is 

SceT 1 1t thC “ d ° f ‘ hiS (« have given theouTpufin 

and the 10 ZmlTtfetionT ‘onctTfl “7^' °' “* 15 ■“ 

K and the azimuth angle * can be calculatS as ’ “ d «" 

_ arctan(^l) 

2 

a 2 = cos( 2 </>)/(A ~C) + (A + C) cos( 2 <£) 
b 2 = cos ( 2 <f>)/{C -A) + (A + C) cos(2 <t>). 

U httht 5 S d Uti ° n req “ ire \ much less computation than the full 15 moment calculation 
o bTc , T, T f ' C ' n,r ° id of the ellips<! must be known accurately for the fit’ 
“7 r n th,S S °‘“ ti0n WiU n °‘ ” k if lhq points that are to be fit represen t ay 
Ihft tLse T b “ aUSe by Summin * the ls ‘ ocder moments the centroid 2 the effipre 
la we points represent is not determined accurately. The threshold/ tracking detector 

t fom a ,ci ° sed ioop ’ - * «-» 2 zt 

the threshold /tr j. dlda ‘ e; lh “ we can be assured that any track that is accepted by 
threshold/tracking routine will be fit well by the 10 moment solution (see appendix B) 
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2 2 -t > 

M22 *M04 *M40 ) *M20 - (Mt3 M22*M04)‘M30 * 2*(MO*M22 M3 t • MO-1 ) • M2 I ♦ M30 - ( M3 I - M2 2 • M V . I * M I 2 • ( M2 2 

2 2 
• M04 *M40 ) *M2 I >*M01 ((MICM3I ♦ M02*M40 - 2*M20*M2? » *M2 1 ‘ <2*M11*M22 - M02‘M3I - M20 • M 1 3 ) • M30 ) * M 1 2 

2 2 2 

• ((MII*M22 - 2*M02*M3t ♦ M20 • M 1 3 ) •M30*M03 » (M13 - M22*M04)*M20 - <MO‘M31 - M22 l*M02*M20 (MI3‘M40 M3 I • 

2 2 2 
M22CMIMM02 (M22 • MdfWOl'MIl )‘M21 - <<M1I*M04 M02*MI3)*M30 - (M02*M40 - M20‘M22 ) ‘M03 ) *M2 1 (( 

2 ? 2 
M I 3 *M3 I M22 ) *M20 (MI3*M-10 M3I-M22MMII - (M3t M2 2 • MIO ) • M02 ) * M03 *M20 <MII*MI3 M02 *M2P > *M30 • 

3 3 2 

MOl • (Mil *M40 - M20*M3 1 ) *M 1 2 (M02*M22 - M20*M04)*M21 )/(2«(((MO M22*M04)*M1I 4 (MI3*M22 M3 1 • M04 ) • Mo? 

2 2 2 
• MD3 *M3 I ) 4 M30 • ( ( M 1 3 *M3 I M22 )*M20 ♦ <MO*M40 M3 1 *M22 ( *M 11 ) *M03 ) ‘M2 I * 2*(((MO‘M3l M22 )*M(12 

2 

- ( M 1 3 *M22 - M3 1 *M04 ) *M I I ) *M30 - ((MO*M3t - 2*M2P 4 M04*M40) 4 M1I - (MO*M22 - M31»M041*M20 * (MI3*M10 

2 

- M3 1 *M22 ) *M02 ) *M2 I ((M0 4 M40 - M3I»M22)*M20 - (M3I M22 *M40) *M I 1 ) *M03 * (M30*M04 - M03‘M31|* 
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M2I • M2 1 *M O » M30 *M03 *M 1 3 ) *M 1 2 - (2*(M30*M13 - M03‘M40)*M21 - 2*<M13*M40 - M3I«M22I*MII • < M3 I M22* 
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Appendix D: Duplicate TVack Elimination 

mOS ‘ di T 1 ’ ° r ‘ brUtC for “’ Way to problem is using an N> algorithm 

th ' * “ d y coordinates of ev^y track record in a file to tho^7JZ 

lone T| d ' ° W<,Ver ’ for a ffle containing 50,000 tracks, this method takes prohibitively 

frame’s position • Tb T redUCe ^ tUne ’ S *° take advanta S e of the sequentiality of each 
T a T TV? he P ° wer scan fiIe - If lhe number of rows and columns scanned is 

"a7tord o ’ T* " lhe frame ind “ n,,,nber “ induded “ “ element in the 

of tracks in ^^.74^77^ C °° rdin8teS ° f ° n,y *° tho “ 

,7° “ sure lha ‘ ^ necessary frames are checked for duplicates, we search the current 
frame the subsequent frame to the right or left (depending on which way the scan pro 

S'the thr f 81 "* 8 n °‘ T‘ St if the CUrrent frame marks the end of a row), and each 

cs, F77“, mOSt b r ath * he *— < a *™ subject to discretion at 

1,7 l ? h he a 0 " 1 * pattern we use for a frame in a row with an original 
search dmecfon from left to right. In order for this scheme to work, tracks TusTbe ta 

“7 : h 't 7° WS d ‘ reCt ^ *o frame number. If no duplicates are Zd for a 

fo a »d r 7itZ 7777 new , ffle ; if a dupbcate is found ’ both tracks are &«» written 

anH ^ ^ d hCn 0nIy tHe S6COnd ° f the of the two (sequentially) is retained 

— — — » — - 
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Figure 1. Schematic diagram of the microscope system. 
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Figure 2. Diagram of the menu system and associated tasks 
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Figure 3 . Enhancement /thresholding edge detection sytem 
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Figure 4a. Laplacian Mask 
Figure 4b. Smoothing Mask 
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Figure 6. Track radius histogram for sheet #38 (front) 
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Figure 9. Plot of Ri vs. R 2 for sheet #38 (front) and sheet #36 (front) 
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Figure 12. Histogram of track minor axis, balloon data, sheet #182 
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Figure 14 . Histogram of zenith angle, balloon data, sheet #182 
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Figure 16 . Results of Laplacian operation on an ideal edge 
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Figure 18 . Frame definition for duplicate track elimination 
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